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Abstract 

Today’s computational methods enable the determination of forces in complex systems, but 
without field validation data, or feedback, there is a high risk of failure when the design envelope is 
challenged. The data of Childs and Bently and field data reported in NASA Conference Proceedings 
serve as sources of design information for the development of these computational codes. Over time 
all turbomachines degrade and instabilities often develop, requiring responsible, accurate, 
turbomachine diagnostics with proper decisions to prevent failures. Tam et al. (numerical) and 
Bently and Muszynska (analytical) models corroborate and implicate that destabilizing factors are 
related through increases in the fluid-force average circumferential velocity. The stability threshold 
can be controlled by external swirl and swirl brakes and increases in radial fluid film stiffness (e.g., 
hydrostatic and ambient pressures) to enhance rotor stability. Also cited are drum rotor self-excited 
oscillations, where the classic “fix” is to add a split or severed damper ring or cylindrical damper 
drum, and the Benkert-Wachter work that engendered swirl brake concepts. For a smooth-operating, 
reliable, long-lived machine, designers must pay very close attention to sealing dynamics and 
diagnostic methods. Correcting the seals enabled the space shuttle main engine high-pressure fuel 
turbopump (SSME HPFTP) to operate successfully. 

Keywords: Turbomachine, seals, dynamics, antiswirl, CFD, models, fluid dynamic forces, rotor/ 
stator clearances 


Introduction 

Operating turbomachines by nature are unsteady systems. The fluid streamlines are interrupted, 
redirected, compressed, diffused, heated, and expanded by fixed and rotating discrete components 
with unsteady heat-release reactions. These same flows are also a source of rotor vibrations, which 
perturb the core and secondary flows, and in turn, propagate disturbances to other parts of the 
turbomachine. As an example, the NASA Energy Efficient Engine (E 3 ) aeronautical gas turbine 
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(Davis and Stearns, 1985) has 32 blades and 34 stators in the fan, 96 blades and 140 stators in the tenth 
stage (last) of the compressor, and 76 blades and 46 stators in the T1 turbine (first stage of the high- 
pressure turbine (HPT)). All provide very rapid, nearly periodic, high pressure and temperature 
variations in the flow field, and all require interface sealing and cavity load balancing with extreme 
care given to control of vibrations. 

Gas turbines and rocket engine turbomachinery have similar interfacing, but the combination of 
cryogens, rapid startup, high rotor speeds, propellant-lubricated bearings, and dynamics are more 
severe in rocket engines. To this one must add noncentered orbitals and unsteadiness induced through 
eccentric thermomechanical loadings and transients associated with the vehicle flight profile (see 
part 3). Figure 1 illustrates the kinematics of the aeronautical gas turbine seal tip over the acceleration 
and full-power phases along with shutdown and soakback of a simulated flight profile, but without 
flight dynamics (engine on test stand). 

In the past turbomachine rotordynamics have contributed to or been the direct cause of many 
design failures. With the maturation of the gas turbine and the demanding high-power-density rocket 
engine turbomachinery, most machines run stably. The exceptions arise when designers are either 
unaware or forced into a situation where the demands of the design outweigh the risks. In these 
situations the potential for failures becomes imminent whenever the design envelope is challenged 
or expanded. However, over time even the best of designs degrade with use, and many turbomachines 
in the field have encountered instabilities due to their operational time and the limited life of their 
components (see appendix B). As a result, it is generally accepted that no two engines operate the 
same or have the same time-variant signatures, even though they exhibit common fundamentals. 
Consequently, responsible, accurate turbomachine diagnostic techniques are required to identify 
root causes of impending malfunctions, with proper operator decisions leading to their prevention 



Figure 1 . — Enlarged view (scale 50:1) of area A. Locus of high-pressure-turbine 
(HPT) blade tip seal relative to static seal during acceleration and full-power 
trip conditions. (Stewart and Brasnett, 1978.) 
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(Bently et al., 2002). However, these important issues will not be addressed herein as D.E. Bendy’ s 
comprehensive book provides very adequate coverage of the subject with many examples, and we 
recommend the reader acquire the book. 

Herein we review the inherent unsteady nature of the turbomachine with small-clearance 
components, specific operational problems of two turbomachine classes and their causes, fixes for 
instabilities, and the predictive methods of Tam et al. and Bently and Muszynska. In part 1 we review 
some sealing requirements and limitations from the viewpoints of the customer, designer, and 
researcher and the controls on the market. In part 3 we examine turbomachine internal flows, engine 
externals, and component life-cycle characteristics. 
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Nomenclature 

amplitude of rotor precession 

normal velocity components induced at segment i by unit vortex strength along 
line segment j 

coefficients in Prandtl mixing-length model 

specific heat 

average swirl speed 

bearing (or seal) radial clearance 

fluid radial damping; turbine diameter 

mean blade diameter 

rotor damping 

transformed space volume element (+1) 
eccentricity 

unit vector 
closing force 

dynamic force due to pressure 
force in seal gap 
spring force 

fluid dynamic radial and tangential forces, respectively 

equivalent radial and tangential forces in Cartesian coordinate system (x, y) 

fluid force 

gravity vector 

blade height 

minimum clearance 

scale factors in three coordinate directions 
determinant of Jacobian matrix 
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fluid radial stiffness coefficient 
rotor or fluid direct dynamic stiffness 
rotor modal stiffness 
external stiffness 

rotor or fluid quadrature dynamic stiffness; lateral spring force coefficient 
fluid tangential stiffness coefficient (cross-coupled stiffness), F f !z 
fluid film stiffness 

shaft partial stiffnesses, left and right sides, respectively 

friction coefficients on rotor and eddy chamber or labyrinth chamber sur- 
faces, respectively 

turbulence energy production-dissipation 

blade length 

seal or bearing length 

mixing-length scale of seal 

rotor modal mass 

fluid inertia 

disk and journal masses 
mass flow 

number of labyrinth chambers 

mass of unbalance 

number of distributed line segments 

unit normal vector 

rotating force 

static pressure 

inlet and outlet axial pressures 
axial pressure drop, P- m - P QUt 
fluid pressure 

source flow; carrying force; lateral force 
nondimensional force 
dynamic pressure 

seal rotor or bearing radius; Richardson flux number 
shaft radius 

fluid local radius in bearing or seal clearance; radius of imbalance 
displacement vector 


NASA/TM— 2004-21 1991/P ART2 


4 


s 

T 

T 


1 total 

t 

u 

^abs 

u 

u 

u,v 


V 

Ko fi 
Volcell 
x, y,z 

Y o 


z r 

0 

a 

P 

y 

8 

8 

8 

C 

11 

e 

r 

x 

X e 


source terms 
torque 

inlet temperature 

torque at average tip clearance 

torque at maximum tip clearance 

total torque of all turbine blades 

time 

fluid relative velocity vector 
fluid absolute velocity vector 

reference frame velocity vector 
local velocity 

reference frame circumferential and radial velocity components, respectively 
preswirl circumferential velocity at radius r 
velocity at midpoint of segment i 
physical volume of grid cell 

shaft horizontal, vertical, and radial complex displacements 
rotor offset 

axial Cartesian coordinate in bearing or seal; shaft radial complex 
displacement, x + jy 

radial displacement magnitude of z 

null matrix 

phase; (T^T^/Tq 

destabilizing stiffness parameter, a/(28 RF f ) 

ratio of specific heats; ratio of hole surface area to seal surface area 
local vortex strength 

average damping coefficient; displacement 
eccentricity ratio, e/c 
contravariant base vector 

normal distance between housing and shaft (stator and rotor) surfaces 
purge flow efficiency parameter; viscosity 
angular position of tangential force 
vortex strength 

fluid average circumferential velocity ratio 
preswirl-related fluid average circumferential velocity ratio 
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Superscripts: 


external fluid average velocity ratio 

injection-related fluid average circumferential velocity ratio 
nonlinear function of fluid average circumferential velocity ratio 
turbulent viscosity; rotor orbital velocity 
fluid kinematic viscosity 

transformed- space arc lengths in three coordinate directions 
fluid density 
stress tensor 

angular coordinate in bearing or seal 

rate of strain in Newtonian fluid; angle measured from position of minimum 
clearance 

load factor 

nonlinear stiffness and damping functions 

rotation; onset speed 

rotor speed 

critical speed 

whirl 

rotor precession (perturbation) frequency, or angular velocity 
rotational speed 

rotor stability threshold 

horizontal and vertical oscillation, respectively 
gradient 


denotes vector 
denotes tensor 


Nomenclature for Thomas (fig. 4): 

A flow area 

C* ratio of flow kinetic energy to shroud pressure difference, 

PQ \I2Apb 

C lu tangential component of absolute outlet velocity of stator blade flow 

A h is stage isentropic enthalpy difference 

K 2 clearance excitation coefficient, ql/U is 

l unit length of rotor blade 
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stage mass flow 

pressure difference at the shrouding 

clearance excitation cross-force spring coefficient, q s + q D 

velocity 

isentropic tangential force, mAh-Ju 
average circumferential velocity 
eccentricity 

pressure distribution portion of q, ^ J sin (p dA 
clearance loss portion of q, cos cp dU 
density 

circumferential angle 
load factor, 2A h is /u 2 


Nomenclature for Benckert (figs. 8 to 10 and eq. (1)): 
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Po 


entry swirl speed 
eccentricity 

ratio of cavity to rotor height (cell depth) 

tangential spring force coefficient 

seal length 

number of chambers 

ambient pressure 

torque or tangential force 

upstream pressure 

lateral force 

rotor radius 

cell width 

rotor speed 

relative eccentricity, e/A r 
upstream density 


Glossary: 

We include a glossary because various authors use the same name for different phenomena and/ 
or different names for the same phenomenon. 


Whirl Whirl is used in general as a synonym of rotor precession or (better) lateral orbiting. Whirl 

is also used to describe specific self-excited vibrations of the rotor, such as fluid whirl. 
Also oil whip, steam whip, aerodynamic whip, fluid-induced instability, clearance 
excitation, half-speed whirl, and probably even more names, describe the same phenom- 
enon: self-excited vibration of rotors due to forces in a fluid environment. 
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Swirl 


Swirl in seals represents circumferential flow in the direction of rotation. The flow is, of 
course, three dimensional and might be unsteady, but the existence of swirl means that in 
the rotor/seal clearance a significant circumferential flow has been established. Active 
antiswirl injections, as well as passive preswirl (constant tangential vanes at the entrance 
to the seal) techniques, clearly were developed to prevent this circumferential flow in the 
seal from being significant. 

X The product of fluid circumferential velocity ratio times rotational speed X(ti R is the speed 

at which the fluid force rotates inside the clearance (i.e., the ratio at which the fluid force 
rotates relative to the rotor speed) and is especially important in the damping force because 
the fluid inertia force may rotate at a different rate. If the circumferential flow in the 
clearance is modified by an external source requiring application of some outside energy, 
X may become larger than Vi, or even negative. For example, if there are injections (such 
as antiswirl in the direction against rotation), the circumferential flow is modified. If this 
external source is very strong, to the point where the circumferential flow is in the 
direction opposite to rotation, X may become negative; however, it is then probably 
inconvenient to relate it to the rotational speed. If the injections are in the direction of 
rotation, X becomes greater than Vi, depending on the amount of injected fluid. Note that 
X higher than Vi occurs normally in pumps with high fluid recirculation. In such cases the 
self-excited fluid whirl has frequency higher than V 2 rotor speed. 


Cases of Small Clearances Between Rotating and 
Stationary Parts 

Seals and Bearings 

Although the primary function of a seal is to control leakage, a secondary but equally important 
purpose is to provide (or at least not to infringe on) rotordynamic stability. Bearings also have a dual 
role, to support the rotor load and to provide dynamic stability. Seal and bearing combinations can 
be used to enhance stability as well as to decrease bearing loads (von Pragenau, 1992; Munson et al., 
2002; and Bently, 2001). Von Pragenau’ s work was so good that Rocketdyne patented the concepts. 
More recent work on fluid film bearings was done by Dimofte and Hendricks (2001), Bently (2001), 
and Bently et al. (2002). 

Further, the flow fields in bearings and seals have similarities, although with some distinct 
differences. The axial pressure drop in the seal is large, and the circumferential pressure drop is 
nominally smaller with coupling to structure, cavity flows, and the power stream (see part 3). For the 
bearing the opposite nearly holds true except that coupling to the power-stream flow is indirect. 
Further, zones of secondary flow appear inevitable. Cavitation phenomena can occur in both 
bearings and seals in the form of either gaseous cavitation (pressure lower than atmospheric) or 
vaporous cavitation (pressure lower than saturation pressure). For these complex flows local bulk 
flow model physics are inadequate, but in many cases the model adequately predicts system 
dynamics and has been successfully used by many authors (e.g., Muszynska, 2001 ; and San Andres 
and Childs, 1997). Seal or bearing design using computational fluid dynamics and knowledge-based 
methodologies can enhance turbomachine stability. (See part 1 and computational methods and CFD 
modeling in part 3). 

Simplifications that incorporate first-order physics are often used in modeling complex flows. 
In this respect a more complete modeling of the turbomachine operating conditions, loads, and 
supports is sought because rotor (component) resonances can occur far from critical and single-point 
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sensing can fail to warn of impending failure (Queitzsch and Fleming, 2001). Although their 
examples are related to rolling-element bearings, these same principles are applicable to interactions 
between small-clearance fluid films and system dynamics throughout the turbomachine. 


Seal/Structure Interaction 

When the clearance gap varies with time, lateral forces occur that act out of phase with case 
distortion. Depending on changes in gap in the flow direction, excitation or damping of the lateral 
forces occurs. The rotating drums shown in figures 2 and 3 represent very large energy storage 
systems that must be properly contained. For example, a 356-kN- (80-klb-) thrust engine stores 8 MJ 




Figure 2. — Rotating drums, (a) Generic gas turbomachine showing low/high dual-spool drum rotors 
and support, (b) Energy Efficient Engine, showing high-pressure turbine. (Halila et al., 1982.) 
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Clearance excitation coefficient, K 2 ; load factor, iJj 



Figure 3. — Steam turbine stage types investigated. (Leie and 
Thomas, 1 980.) 
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designs (A, B, C) 
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A — Integration of pressure 
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Calculation from 

leakage losses 


(b) 


Figure 4. — Steam turbine instabilities, (a) Three labyrinth seal configurations investigated, 
(b) Test results. (Leie and Thomas, 1980.) 
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Figure 5. — Compressor discharge annular 
diffuser. Where compressor discharge seal 
is at a substantially smaller diameter than 
inner wall of main flow annulus, check for 
possible cross flow and cellular recircu- 
lating flows between wheel space chamber 
and diffuser inlet. (Alford, 1967.) 



n = 1 circumferential waves n = 3 circumferential waves 

Figure 6. — Wheel space chamber coupled to diffuser 
inlet. Where an axial gap connects wheel space 
chambers and diffuser inlet, cellular recirculating 
flows may occur. (Alford, 1967.) 


in the high-pressure turbine (HPT) rotor system and 18 MJ in the total engine; a 445-kN- (100-klb-) 
thrust engine stores 18 MJ in the HPT rotor system and 60 MJ in the total engine. 

Early on, Smith (1933) studied the stability of rotating systems and reported that stationary 
damping in bearing supports always favors stability but that damping in rotating parts should be 
avoided when the machine has to run above a critical speed. He believed that unsymmetrical flexible 
bearing supports may improve running (asymmetry is to be avoided in shafts and rotors) and 
cautioned not to run at speeds close to critical or to half one of the normal frequencies of the system. 
Sometimes this anisotropy is unavoidable (as in two-pole generators). 

Some 25 years later self-excited vibrations engendered by forces of the fluids being processed 
were investigated by Thomas (1958) (see also NASA CP-2133, 1980, p. 303, fig. 3), who looked into 
steam turbine whirl excitation. For a bearing, in certain respects, stiffness and damping in themselves 
are not so important, but their relation to one another is important. Clearance excitation forces emerge 
from nonuniform flows (leakages) and pressure distributions as a function of eccentricity. Thomas 
introduced swirl-preventing sheets of metal that reduced or eliminated the excitation forces (fig. 4), 
where the clearance excitation coefficient K 2 and load factor \[ / are plotted against the ratio of flow 
kinetic energy to shroud pressure difference C E . 

Alford (1963, 1965, 1967) followed with a series of papers establishing some methods of 
controlling instabilities in aeronautical turbomachines. All or nearly all unstable and oscillatory 
flows are self-excited by flow systems with rising pressure characteristics (one in which the internal 
pressure increases with an increase in flow), such as every cycle in an eccentric seal, housing, 
bearing, or compressor (fluid-handling machines in general). Tangential air velocity has a primary 
effect on resonances with circumferential waves. Phase differences between longitudinal waves in 
parallel diffuser passages give rise to spinning pressure waves in adjacent annular ducts (figs. 5 
and 6). Resonance occurs when the space oscillations and spinning waves are the same for the wheel 
and disk (reminds one of Ben Franklin’ s glass harmonica). Alford employed principles of designing 
high-frequency waveguides in his analytical work. 
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Stator 



Figure 7. — Labyrinth seal vibrations. (Armstrong, E.K., discussion 
of Alford, 1967.) 


Self-Excited Seal Drums 

Description . — In a discussion of Alford’s paper (1967), E.K. Armstrong cited cracking in the 
compressor-discharge-pressure (CDP) labyrinth seal stator drum. In this case the stator vibrated to 
1300 Hz (perhaps similar to a more recent engine failure problem) with 6 nodal diameters and 
12 nodes with a distortion (fig. 7). At a certain temperature the frequency of mode vibration in the 
annular air cavity with six wavelengths becomes equal to the sixth diametrical mode of the stator 
(7TD/6). Under these conditions variation in the air gap in the labyrinth due to the stator vibration 
excites vibration in the cavity. In turn the resulting pressures act upon the stator, and if the phase is 
correct, self-excited oscillations occur. The mechanism could be stopped by inserting 24 equally 
spaced radial baffles across the air cavity (phase and baffles are cited by Thomas, 1958). Alford 
comments that the mechanism that excites the labyrinth seal stator appears to be the same 
phenomenon that excites the labyrinth seal rotors. 

Below are “rules of thumb” that exist in the field. These rules are listed here without specific 
classification and without any evaluation for applicability. Therefore, they are often misused. 

Some rules of thumb . — Prevent parasitic leakages. Provide radial baffles near the wheel rim to 
prevent cellular recirculation. Provide a cylindrical surface adjacent to the baffle to ensure that a 
minimum area of the flow communication slot is independent of axial motion. Reduce swirl entering 
the diffuser. Provide a margin greater than 15 percent between resonant longitudinal waves and 
wheel space waves. If inlet flow is split, extend the duct splitter to the compressor inlet to prevent 
crossflow between inlets. Provide adequate radii on passage leading edges to prevent vortices and 
flow disruption, particularly flow separations. Check rotor blade excitation modes for excitation and 
avoid them. Dynamically monitor parameters (e.g., with proximity probes). 

The work of Benckert and Wachter (1978, 1980) set the stage for further understanding the basic 
forces in turbomachines and the advancement of seal flow inlet swirl as a major source of high- 
power-density turbomachine seal instabilities. 


Swirl Brakes 

Although known and reported by Thomas and Armstrong, from the work of Benckert and 
Wachter came the widely acclaimed swirl brake, which did solve several labyrinth seal instability 
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Recessed casing 


Mortised rotor: 
m = 10, 6, 4, 2 
A r = 0.25 mm 


(c) 



Mortised seal strips 
in casing and rotor 


m = 23, 11,7, 3 
A r = 0.5 mm 


Testing parameters of labyrinth: 

Entry swirl speed, C u q Rotor speed, U w 

Differential pressure, Ap s t Number of whirling 
Eccentricity, e chambers, m 

Geometry, Ar, h, t 


Figure 8. — Close-clearance labyrinth shaft seals, circumfer- 
ential types. (Benckert and Wachter, 1978.) 


problems in U.S. industrial power turbines. In their work Benckert and Wachter (1978) studied 
( 1 ) teeth on stator, smooth rotor; (2) teeth on rotor, recessed stator; and (3) interlocking teeth on rotor, 
teeth on stator (fig. 8). They found that lateral forces related to circumferential flows are detrimental 
to rotor stability and that blocking the flow could control these perturbations (fig. 9). Here the lateral 
spring force coefficient Kq is given as a function of rotor relative eccentricity 8; Kq can be estimated 
from the swirl speed and for see-through labyrinths 


2 (K e «r 2 ) 2 « p 0 cj; 0 (P 0 ~ p a> 


( 1 ) 


For more established flows 

\K Q (htr)M = - k 2 ) c 2 „+ (2) 

where the subscript m is the average value in the labyrinth chamber, m n is the number of labyrinth 
chambers, k 1 is the friction coefficient on the rotor, k 2 is the friction coefficient on the eddy chamber 
or labyrinth chamber surfaces, and C um the average swirl speed. (See also Muszynska’s lumped- 
parameter model in part 2 of appendix A.) 
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Figure 9. — Speed -induced lateral forces by using swirl webs in each 
eddy chamber of labyrinth. Tangential spring force coefficient, 
P a /Po = 0.49, h = 2.75 mm. (Benckert and Wachter, 1980.) 


According to Benckert and Wachter (1978) the swirl brake work was an outgrowth of a 
description in Wohlrab (1975, translated in 1983) where a steam turbine was stabilized by placing 
the brake upstream of the first tooth to redistribute the entry flow in the seal and control the 
circumferential flow in the chambers. Benckert and Wachter also studied two sizes of honeycomb, 
4.95-mm (0.195-in.) cell by 3.75-mm (0.148-in.) depth and 3.55-mm (0.140-in.) cell by 2.75-mm 
(0.108-in.) depth. They showed relative lateral (destabilizing) forces an order of magnitude larger 
than for the interlocking labyrinth (fig. 10). Iwatsubo and Iwasaki (2002) evaluated a swirl brake and 
a labyrinth seal similar to those of Benckert and Wachter. Iwatsubo and Iwasaki found that adding 
control fins in the labyrinth grooves was effective in reducing the circumferential swirl velocity and 
that adding control fins at the inlet, particularly the first and second grooves, controlled seal stability. 
This effect was predicted by Tam et al. (1988) (see the section Computational Method of Tam et al. 
and part 1 of appendix A). In yet another method termed a web seal (appendix B), slots are introduced 
to enhance stability. 
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Figure 10. — Lateral force spring coefficients as function of relative 
rotation flow value for two honeycomb seals. Seal clearance, 
r, 150 mm; A r, 0.25 mm; seal or bearing length, /, 99 mm. 
(Benckert and Wachter, 1 980.) 


Pumps and Impellers 

Henry Black contributed heavily to the rotordynamics field. Black (1968) analyzed annular 
clearance spaces of impeller and balance piston seals. He found large damping associated with 
squeeze of the fluid to be comparable to stiffness. With dominant inertia effects, apparent negative 
stiffness can occur (see appendix A). Analysis gives synchronous critical speeds, whirl amplitudes, 
and whip boundaries. Fluid rotation in synchronous whirl reduces the damping forces. At high speeds 
fluid rotation may lead to whip, which can be avoided by tightening the clearances but could occur 
as the seal wears. 

Thompson (1980) analyzed 115 individual compressor stages contained in 20 multistage bodies 
and concluded that the average damping coefficient 8 per stage to ensure stability is 

^total = 1.85 = S fluid + S mechanical and S fluid = -F x !u t (3) 

where F x is the cross-coupled force component (tangential force) due to dynamic excitation and u t 
is the rotor orbital velocity. Thompson also cites the opportunity to add damping to avoid the effects 
of compressor surge. 

Colding- Jorgensen (1980) modeled the centrifugal compressor or pump stage with a vaneless 
volute based on a two-dimensional representation of the diffuser and impeller by an equivalent vortex 
of strength T and a source flow Q at a point. The velocity at any point is the sum of the velocities 
induced by singularities in the field 
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lAijV = v -j"i 


(4) 


where • is the velocity at the midpoint of segment i, rif is the unit normal vector, and A- are the 
normal velocity components induced at segment i by unit vortex strength along line segment j. The 
sum is over the N distributed linear segments along unit length line segments j covering the volute. 
The N sources at the midpoints are the unknown local vortices of strength y induced at point i. The 
velocity V ^ • is induced by the vortex (T, Q ) representing the impeller at a given eccentricity and 
velocity of the vortex source. Once the solution for Jj is achieved, the velocity at any point in the 
plane can be found. (See also NASA CP-3344, p. 177.) 


Effect of Damping 

Although it was noted by Smith (1933) sometime earlier, Crandall (1980) provides a limited 
explanation of how rotating damping can destabilize a rotating system. When rotation Q is faster than 
the whirl co^, rotating damping drags the rotor along; and when slower, it retards the rotor. When 
stationary damping is present, the stability line is balanced by backward and forward drag forces. At 
low rotational speeds the rotating (internal) damping force in rotating elements acts collinearly with 
the external damping force (damping in nonrotating elements and the rotor environment); at high 
rotational speeds the rotating damping force changes the polarity — thus, it opposes external damping 
up to its nullification. A forward whirl at rate co w with respect to stationary axes appears as backward 
whirl at a rate of Q - co^ with respect to a system rotating supercritically at rate Q. Growth or decay 
is estimated from the energy balance. If the system is horizontal and elastically supported (e.g., 
compliance is greater in the x direction than in the y direction), the system becomes a superposition 
of a horizontal oscillation co 1 and a vertical oscillation co 2 >0)^ The more anisotropic the supports, 
the greater the difference between co 1 and co 2 , an d there is no tendency to set up a whirl of the type 
that can be dragged forward by rotary damping until these forces have been increased by risings in 
speed to be commensurate with the difference in elastic forces in the two principal directions. Thus, 
anisotropic bearing supports should raise the onset speed Q for instability due to rotating damping. 
Nevertheless, the pendulum in the rotating coffee cup explanation for a fluid film bearing when its 
rotational speed exceeds 2x (twice the critical speed of the rotor) is also interesting. 

As a possible illustration Crandall considered a thin fluid layer in a cup rotating at Q, a pendulum 
in the fluid layer rotating at 0d n (both in the counterclockwise direction (fig. 11)). When Q< co n , fluid 
drag retards the pendulum bob and the radius of the orbit decreases. On the other hand when Q > co w , 
fluid drag pulls the pendulum bob around in an orbit and the radius increases — a possible explanation 
for whip when Q > 2co^. Here, Crandall assumes that X = 0.5. The more general assumption that X 
is not constant and, in particular, differs from 0.5 opened new insights into the fluid forces acting on 
rotors in relatively small clearances (see part 2 of appendix A). Note that for the same rotor model 
the fluid whip condition used in the Bently-Muszynska fluid model is Q = co^ IX. Because X can be 
modified (e.g., by swirl brake, fluid injection, or applied pressure), Crandall’s instability threshold 
also can be modified. An example is given in part 2 of appendix A. 

For a more complete assessment of theory, data, and analysis in the development of seals-related 
instability problems and their solutions and seals and secondary flow codes development, see a 
NASA Conference Publication as given in the references. For further work in the area of pumps and 
turbomachines contact 

Dara Childs, Texas A&M Pump Consortium Meetings, College Station, Texas 

Allen Acosta and Chris Brennen, CalTech, Pasadena, California 

Yoshinobu Tsujimoto, Osaka University, Osaka, Japan/ ISROMAC-9, Hawaii, 2002 
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Figure 1 1 . — Rotating coffee cup 
illustration. (Crandall, 1980.) 


Bently Pressurized Bearings Company 
NASA Conference Proceedings 

It was from these and other field experiences that the space shuttle main engine turbomachinery 
field came to realize the importance of lateral forces that destabilize the rotor and the importance of 
upstream correction of the inlet swirl. In some seals the swirl brake mitigates those forces. The effort 
to control circumferential flow continues to this day (see appendix B). 


Thrbomachines With Specific Operational Problems 

High-Power-Density Turbomachines 

Let us digress for a moment as some background seems in order. For more complete work, see 
books by Childs (1993) and Vance (1988); and for diagnostics see Bendy et al. (2002). As noted, 
steam turbine seal leakages were investigated by Thomas and later by Alford for aeronautical turbine 
rim seals. Although the destabilizing forces of bearings were somewhat recognized and resolved, the 
forces in seals of the space shuttle main engine high-pressure fuel turbopump (SSME HPFTP) were 
not. Instabilities in the SSME HPFTP caused many failures and prompted Otto Goetz (chief 
engineer) and his group (e.g., George von Pragenau, Henry Black, Ed Gunter, Larry Ludwig, Dara 


NASA/TM— 2004-21 1991/P ART2 


17 



Figure 12. — High-pressure fuel turbopump (HPFTP) 
shuttle seal, (a) Three-step labyrinth, (b) Straight 
cylindrical, (c) Three-step cylindrical. (Hendricks, 
1981 a,b,c.) 



X 


Figure 13. — Damping seals and bearings. 

(a) Damping seal concept with isogrid 
stator, (b) Damping bearing concept with 
stationary orifices (orifices can also be in the 
shaft). (Inventions of George von Pragenau, 
Huntsville, Alabama; Hendricks et al., 1995.) 
See also part 1 , figure 33. 



Childs, and others) to propose investigating the seals as a probable cause. Goetz suggested full-scale 
testing of the HPFTP seals in a simulated, nonrotating configuration. Hendricks tested and later 
published dynamic force data for three seal configurations (labyrinth, cylindrical, and stepped 
cylindrical) with cryogenic and gaseous hydrogen as the working fluid for a range of conditions. 
These data along with the stability and damper seal concepts of von Pragenau (1982, 1985) served 
as guidelines for subsequent turbopump refurbishment (Hendricks, 1981a,b,c) and kept the space 
shuttle “flying” until a better understanding emerged (figs. 12 and 13). Correcting the seals enabled 
the SSME HPFTP to operate successfully. 
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Figure 14. — Hole patterns for damper seals, where 7 is ratio of hole surface area to seal 
surface area. (Yu and Childs, 1997.) 


These turbopump instability problems led Prof. Childs to organize the first seals and rotordynamics 
workshop at Texas A&M as reported in NASA CP-2 1 3 3 (in retrospect, a classic document) . A dearth 
of data and a wealth of seal instability problems coupled with little understanding of sealing forces 
in high-power-density turbomachines dominated the workshop. As understanding began to emerge, 
the results were documented in NASA Conference Publications (see references). A sampling of 
topics include field test data, seal dynamic coefficient measurements, analytical methods, numerical 
techniques, and practical applications. Computational fluid dynamics (CFD) has merged the 
classical dynamics approach with that of the actual flow field, with some success in describing these 
forces. Damping concepts of von Pragenau have been extended to a variety of hole patterns (fig. 14). 
Swirl brakes and antiswirl methodology have been applied with significant success (fig. 15). Design 
data, from foreign and domestic sources, provided dynamic coefficients; and the data of Braun et al. 
(1987) quantified the flow fields near the minimum clearance for both the convergent and divergent 
zones of a simulated bearing. Seals began to play a prominent role in stabilizing turbomachines, and 
some understanding of sealing forces in high-power-density machines began to emerge, although the 
details continue to be analyzed yet today. 

The Benckert and Wachter explanations led rotordynamicists such as Crandall to also give 
explanations of the instability problems; Black, Fleming, Vance, Childs, Gunter, and others provided 
theoretical analyses. NASA initiated two programs, one with Childs dealing with seal rotordynamic 
instabilities and the other to design seals and secondary flow systems. The rotordynamics program 
with Childs provided data, analysis, and biannual workshops that over the years enabled a 
considerable degree of confidence to emerge within the turbomachinery community. NASA initiated 
the development of the computer codes SCISEAL and INDSEAL with progress of the seals 
community reported at workshops on an annual basis. Later San Andres developed TAMSEAL. 
These codes were an outgrowth of the need to capture experience, data, and analysis into code form 
for seal and secondary flow designers. 
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Figure 15. — High-pressure oxygen turbopump (HPOTP). (a) 
Interstage seal (relative dimensions not to scale), (b) Radial 
swirl brake, (c) Improved swirl brake. (Childs et al., 1991 .) 


A long series of experimental seal testing by Childs, Vance, Iwatsubo, San Andres, Bendy, 
Muszynska, and Allaire providing dynamic coefficients, analytical methods (Fleming, 1 Gunter, 
Nelson, Black, Brown, Muszynska, Bendy, and others), and computational techniques (Tam, 
Przekwas, Rhode, Nordmann, Athavale, San Andres, and others) began to replace the data and 
analysis dearth with credibility and served as the basis for codes developments and analysis. The 
CFD numerical work of Tam et al., to be discussed later, in terms of dynamic stiffness, and the works 
of Nordmann that centered on flow CFD and the small-perturbation method with cross-coupled 
terms paved the way for use of CFD modeling with full-field-coupled conjugate solutions developed 
by Athavale (see part 3). 

Spurred by the many industrial field instability problems and the seals data, Bently et al. (1986) 
evolved their work on dynamic stiffness (direct and quadrature) relating the destabilizing circum- 
ferential flows by the swirl ratio X. Through this parameter, which relates the magnitude of the 
circumferential flows along with the rotor forces, phase relations were determined and criteria for 
self-excited vibrations established. The relations of whirl (self-excitation frequency linear with rotor 
speed) and whip (self-excitation frequency independent of rotor speed) in seals and bearings took 
shape (fig. 16). Muszynska (2001) gave a succinct summary of these results in an invited lecture at 


According to David Fleming of NASA Glenn, he and Gunter looked into the HPFTP instability problem and found 
the analysis had been carried out for a cylindrical seal. Had this seal been used rather than the labyrinth seal, sufficient 
damping might have been available to avoid the instability. 
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Figure 15. — Concluded, (b) Radial swirl brake, (c) Improved swirl brake. Dimensions 
in millimeters. (Childs et al., 1991.) 


ISCORMA-1 . In each case the emphasis is related to controlling the lateral or tangential forces that 
destabilize the rotor while maintaining those dynamic (direct) forces necessary to support turbomachine 
loadings (see also Bently et al., 2002). 

In a recent effort Kawak et al. (2002) investigated a seal concept for an SSME-class liquid 
oxygen-methane pump. Held by a clamping nut and inserted between the rotor and the pump case, 
a floating ring seal comes into contact with the rotor at its maximum eccentricity position (fig. 17). 
As the pump comes up to speed, the backflow from the impeller forces the floating ring against the 
static ring while fluid dynamic forces position it in orbit about the shaft center. At this point, the 
floating ring behaves as a cylindrical seal. Their dynamic analysis shows that a floating ring could 
provide pump stability to 80 000 rpm whereas a rigid cylindrical seal could provide pump stability 
to only 40 000 rpm, the latter being unstable in the operating range of 50 000 rpm. 
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Figure 16. — Fluid whirl and fluid whip modes. 
(Orbits filtered to whirl or whip frequency.) 

(a) Whirl mode. (Disk lags journal very 
slightly. Shaft behaves like a rigid body.) 

(b) Whip mode. (Shaft deflection is similar 
to first balance resonance mode but is 
corkscrewed. Disk lags journal by up to 
90°.) (Muszynska, 2001 .) 


Clamping nut 


Supporting ring 




Figure 17. — SSME-class liquid oxygen-methane floating ring pump seal. (Kawak et al., 2002.) 


Aeronautical and Aeronautical-Derivative Turbomachines 

Returning to the perhaps more familiar aeronautical gas turbine industry, Abbott (1981) provides 
a good review of aeroelastic instabilities amid sets of concerns. None of the results can be reproduced 
from these reviewed papers starting with Alford’s work in the late 1950’s and early 1960’s (Alford, 
1963, 1971) and all the subsequent ones. To be fair, the problems and solutions are reported but key 
geometric information is always missing. Each company must have its own files or know the engine 
parameters for those results that are held proprietarily (i.e., this is a guide; get your own test data). 
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failed; stators supported 
at exit have not failed 

For Pb > P A Rotating components 
supported at exit have 
not failed; stators sup- 
ported at inlet have failed 


Figure 18. — Inner and outer labyrinth air seals, (a) Damper ring, (b) Damper drum (Abbott, 1981). 
(c) Effect of seal component support (Ludwig, 1978). 


Damper Rings and Drums 

Engine companies preferring to support seals at the low-pressure side have no reported seal 
vibration problems with this configuration (e.g., outer seal (fig. 18)). However, supporting a seal at 
the high-pressure side highly stresses the “knee” or last tooth region, and instabilities have caused 
engine failures (e.g., inner seal (fig. 1 8)). Here, the classic “fix” is to add a damper ring or a cylindrical 
damper drum as in the case where the ring has insufficient damping. These are usually severed rings 
(also called split rings) or severed drums (fig. 18). Abbott’s (1981) analysis included the circumfer- 
ential fluid velocity that enabled him to show that damping of a seal supported on the low-pressure 
side (outer seal) crosses over from apparent negative to positive damping before the mechanical 
system’s natural frequency reaches the fixed acoustic frequency (fig. 19(a)). On the other hand for 
a seal supported on the high-pressure side the damping becomes apparently negative when the 
mechanical system’s natural frequency exceeds the fixed acoustic frequency (fig. 19(b)). Figure 19 
also gives experimental stability data for low- and high-pressure-side supports. Three different seal 
configurations (no details) were tested to confirm previous engine concerns, but as a safety factor 
they added a seal drum to the seal ring. 

In his conclusions Abbott (1981) states “results are for stationary seals which are designed to 
have either high- or low-pressure-side support. For a rotating seal the mechanical natural frequency 
is the frequency which would be measured by an observer moving at the same tangential velocity as 
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Figure 19. — Seal support and damping characteristics, (a) Aerodynamic damping characteristics, 
(b) Correlation of seal support test results. (Abbott, 1 981 .) 


the air in the seal tooth cavities.” For each mode consideration must be given to both forward and 
backward mechanical wave frequencies. It is clear that circumferential flows can destabilize the seal, 
and according to Abbott’s definition of mechanical frequency it is related to the inverse of the 
^-parameter of Bently et al. ( 1986), where X is the ratio of average fluid-force rotational speed to rotor 
speed. 

Rule of thumb . — With reference to figure 1 8, if a rotating seal is supported on the high-pressure 
side, add dampers. If it is supported on the low-pressure side, there is generally not a problem. Make 
the seal stiff enough to ensure that the mechanical frequency always exceeds the acoustic frequency. 
For engine companies that support their seals at each end, these criteria do not hold. For example, 
Lewis et al. (1978) investigated high-cycle-fatigue cracking in an F-100 engine CDP seal where the 
rotating seal is supported at each end and the static seal is supported at the midspan (fig. 20; similar 
to fig. 10 in part 1). The cracking was common to seals with tight clearances. In rig testing a cavity 
bleed was used to perturb and excite the seal (fig. 21). Engine thrust and specific fuel consumption 
(SFC) were very sensitive to flow through this seal as is usually the case for the CDP seal (Hendricks 
et al., 1994). Stability was 13 times more sensitive to seal tooth 4 (downstream) than to seal tooth 1. 
Although the actual clearances are not given, the fix was to convert it into a converging seal. Opening 
the first tooth clearance by 0.254 mm (0.010 in.) (radial) and closing the fourth tooth clearance by 
0.0381 mm (0.0015 in.) (radial) stabilized this seal configuration, but the actual clearance is not 
given. For a smooth cylindrical seal optimum stiffness is achieved for a clearance convergence ratio 
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Figure 20. — Engine cross section, seal support points, and materials 
specifications (Lewis et al., 1 978.) Seal similar to figure 1 0, part 1 . 



by externally operated 
open/close valve 

Figure 21. — Seal land bleed tube testing, 
amplified instability to confirm predicted 
pressure drop effect. (Lewis et al., 1978.) 


of 1.8 (Fleming, 1977). One cannot assume labyrinth seal wear to be uniform and need to evaluate 
the stability of worn seals. 


Aeroelastic and Aeroacoustic Effects 

Schuck and Nordmann (2002) recently investigated aeroelastic instabilities in the range 0.8 to 
2.9 kHz for a nonrotating, two-teeth-on-rotor labyrinth seal. Flow through their apparatus could be 
in either direction, providing seal support either on the high-pressure side (HPS) or from the low- 
pressure side (LPS). Both the stator and the rotor cylinders were supported at the same end, a 
difference from the configuration of figure 16 (Abbott, 1981), making direct comparisons more 
difficult. For the LPS system the measured cavity pressure amplitudes diminished. For the HPS 
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system the growth of the pressure amplitude was reproducible and higher in the cavity than in the 
labyrinth. In their analysis Schuck and Nordmann (2002) combined the acoustic equations with finite 
element methods to provide solutions for the rotor and stator frequency and mode shapes. Mode 0 
happens but is not as likely. Mode 2 and mode 3 are computed with and without preswirl elements. 
Mode 3 (fig. 22(a)) is difficult to distinguish from mode 6 unless phase is considered. However, 
mode 2 is more readily distinguished (fig. 22(b)). This experimental and numerical work begins to 
lend further understanding to Alford’s rules of thumb and brings potential resolution of the types of 
dampers (rings or drums) and the required geometry into the open literature. 

Schuck and Nordmann (2002) investigated the stiff rotor/stiff stator combination with nominal 
clearance of 0.15 mm (0.006 in.). For high-pressure-side support, mode 2 vibrations appeared to be 
of considerable amplitudes at a pressure ratio (out/in) of 0.86, where inflow pressure was about 7 bar 
(102 psi). For low-pressure-side support, mode 3 was dominant with similar amplitudes, but now the 
pressure ratio to achieve these amplitudes dropped to 0.5 to 0.75 (toward choking conditions) and no 
vibration occurred at inflow pressures of 3 and 5 bar (43.5 and 72.5 psi). Mode 2 vibrations at LPS 
occurred only two times. The HPS and LPS mode 2 results of Schuck and Nordmann (2002, 2003) 
are compared in table I. 

Although there are some differences in measurements, when comparing frequencies one finds 
that LPS is around 690 Hz and HPS is around 890 Hz. There is also a significant difference between 
the pressure amplitudes in the labyrinth chamber (see bold values in table), LPS being nearly half that 
for HPS. The difference between pressure amplitudes in the labyrinth chamber and the cavity is even 
larger for mode 2 than for mode 3 LPS and may be a reason why mode 3 LPS occurs much more often 
and why LPS vibrations need higher inflow pressure and a smaller pressure ratio than HPS vibrations. 

Low-pressure-side-supported structures show less accurate coherence between the structure and 
the pressure field and therefore need higher pressure differences between in- and outflow to achieve 
comparable structural (and pressure) amplitudes. The restriction concerning the applied pressure ratio 
is important. 

According to Schuck and Nordmann (2003) the measured data tabulated in their paper each mark 
an upper or lower boundary for the area of amplitude values (i.e., for a nearly constant inflow 
pressure, all measurements with a certain configuration show amplitude values between those two 
measurements). 


Compressor and Turbine Blade Tip Instabilities 

Under normal operating conditions the leakage flows maximize in the leading-edge region and 
track across the passage to the trailing edge of the preceding blade. The stage-whirl instability is 
instigated (1) when a blade stalls, developing a separation region that unloads the blade and failing 
to produce sufficient pressure to maintain flow around it; (2) when a reduced-flow region forms 
ahead of the blade (blockage), diverting the flow around it and increasing the angle of attack of the 
oncoming blade while reducing it on the preceding blade, and (3) when the differential angles of 
attack cause the stall zones to propagate from the pressure side to the suction side of each blade 
opposite the direction of rotation in pairs or sets of stall zones that can envelop the rotor tip region 
(fig. 23; Johnsen and Bullock, 1965). Blade stall (airfoil flow separation) can occur for a variety of 
reasons, such as surface roughness on blade or casing interface, casing treatment (e.g., felt metal, 
honeycomb, or smooth), tip clearance, tip wear, vortex shedding, foreign object damage, nonuniform 
loading, and eccentric rotation). 

Mailach et al. (2001) cite three types of instability originating from compressor blade tip 
vortices: (1) a local rotating stall region over one or a few blade passages that grows into a rotating 
stall cell, (2) wavelength disturbance that rotates around the annular clearance region and eventually 
breaks down the flow field, and (3) a group of superimposed modes or part-span stall that can 
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Figure 22. — Mode shape of seal housing, (a) Mode 3. (b) Mode 2. 
(Schuck and Nordmann, 2002.) 
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TABLE I.— LOW-PRESSURE- AND HIGH-PRESSURE-SIDE SUPPORT DATA 
WHERE MODE 2 OCCURS 


Support 

side 

Inlet 

pressure, 

bar 

Pressure 

ratio 

(out/in) 

Frequency, 

Hz 

Amplitude, 

pm 

Labyrinth 

chamber 

pressure, 

bar 

Cavity pressure, 
bar 

Rotor 

Stator 

At 

labyrinth 

chamber 

At mid- 
plane 

Near 

support 

LPS 

6.9 

0.75 

684 

45 

10 

0.088 

0.265 

0.25 

0.229 

LPS 

6.9 

0.586 

694 

50 

3.8 

0.07 

0.28 

0.27 

0.248 

HPS 

7.05 

0.9 

881 

37 

52.4 

0.146 

0.25 

0.27 

0.286 

HPS 

7.1 

0.86 

899 

47 

89.1 

0.26 

0.344 

0.372 

0.34 


Direction 



Figure 23. — Propagating stall in cascade. (Fig. 230, 
p. 313, Johnsen and Bullock, 1965.) 


intensify tip clearance noise and blade vibrations. Although not cited, vortex flows within the blade 
passage are engendered by flows along the blade surfaces from the pressure to the suction side, and 
corner vortices may be regarded as yet another class. Under normal operating conditions the leakage 
flows maximize in the leading-edge region and track across the passage to the trailing edge (for more 
detail, see parts 1 and 3). In the experiment of Mailach et al. (2001) third-stage rotor blade data of 
a four-stage compressor were taken at large clearances (4.3 percent ratio of tip clearance to tip chord 
length). These data revealed a strong tip vortex with a reversed flow region within the gap that 
intensified (repeated) every other blade and subsequent blockage. Unshrouded compressor and 
turbine blade tip flows are illustrated in part 1 and blade passage flow computational methods in 
part 3. Hendricks et al. (1995) and Athavale et al. (1993) provide simulated cases. 
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Figure 24. — Compressor map for baseline compressor and compressor A at 
0.0962-cm offset and data set points 3 to 6. (Storace et al., 2001 .) 


TABLE II.— COMPRESSOR A FORCES AND NORMALIZED CROSS-COUPLED 
STIFFNESS (WHIRL ANGLE) FOR ROTOR OFFSET Y 0 OF 0.0965 cm (0.038 in.) 


[Mean blade diameter times blade height, Dfl, 30.3 cm 2 (4.7 in. 2 ); fluid tangential stiffness 

coefficient, K xy = FJY 0 .] 


Compressor 
operating condition 

Test 

point 

Radial force, 
N 

Tangential force 

F y , 

N 

Torque, 

T, 

N-cm 

KJ) p H!T 

Peak efficiency 

3 

-2.71 

2.63 

1010.6 

-1.64 

Design point 

4 

-4.13 

4.36 

1000.6 

-2.76 

High load 

5 

-7.29 

5.69 

989.4 

-4.63 

Peak pressure 

6 

-9.61 

7.88 

926.2 

-6.58 


In discussing compressor whirl Storace et al. (2001) claim Alford (1965) concluded that blade 
loading would be lower in the minimum clearance than in the maximum clearance resulting in 
forward whirl, similar to turbines. Ehrich (1993) concluded the opposite. Citing the work of Thomas 
(1958) and Alford (1965), Storace et al. (2001) statically offset the rotor of a four-stage, low-speed, 
1 .524-m (60-in.), 804-rpm compressor and determined the sum of forces perpendicular to the plane 
of displacement. A sample map is shown as figure 24. 

See table II for a representative set of data forces resulting from the integration of the pressure 
changes about the blade. Figure 25 illustrates these peripheral normalized pressure profiles for the 
high-load case ( P f5 ). 

The resulting forces promoted backward whirl over much of the operating range with islands of 
zero and positive forward whirl near the design point confirming Johnsen and Bullock (1965). For 
the high-load case (P t5 ) the magnitude of the cross-coupled force increased with offset (fig. 26). In 
most cases these instabilities are “fixed” by the manufacturer using split rings, drums, or enhanced 
dampers. However, the root cause remains unexplored and is of little consequence until the design 
parameters call for high stage loadings and marginal compressor designs. 
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Figure 25. — Peripheral static/dynamic pressure ratio P s /q 
variation from leading edge to trailing edge for pressure 
and suction sides of stator blade, (a) At 5% span, (b) At 
50% span, (c) At 90% span. (Storace et al., 2001 .) 
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Figure 26. — Cross-coupled force at high compressor loading Pf 5 . 
(Storace et al., 2001 .) 


Storace et al. (200 1 ) also noted the peaking of the unsteady forces leading the minimum clearance 
by 40° for compressor A. Doubling the clearance at the rotor interface dropped torque efficiency by 
1.6 points versus 0.9 point for doubling the hub stator clearance. The effects of entry cavity are 
explored by, for example, Millsaps in NASA CP-3239 (pp. 179-207). Storace in NASA CP-3122 
(pp. 285-299) discusses predicting the performance of aerodynamically excited turbomachines, but 
Storace et al. (2001) does not discuss the many works appearing in the NASA Rotordynamic or Seals 
Workshop proceedings or in ISCORMA, which would have benefited the authors. 

We will now review and expand the earlier work of Tam et al. (1988) and the lumped-parameter 
model of Muszynska (1986a, b) and Bendy et al. (1986). 


Specifics of Tam et al. and Bently and Muszynska Models 

Computational Method of Tam et al. 

Tam et al. (1988) considered a shaft rotating clockwise at rotational speed top, precessing with 
amplitude A, and having dynamic eccentricity at rotor precession (perturbation) frequency CDp within 
a three-dimensional, static, cylindrical seal or bearing housing (fig. 27). Strong axial pressure drops, 
significant Coriolis and centrifugal forces, and convergent and divergent zones are the major 
contributors to dynamics. For these nonconventional flow fields in very narrow passages with 
passage geometry changing periodically with time, gyroscopic motions are predominant. 

A three-dimensional, rotating coordinate frame attached to, and precessing with, the rotor at 
perturbation frequency C0p can be described. In this coordinate system the given seal and bearing 
geometry is invariant in time (Mullen and Hendricks, 1983; Tam, 1985; and Tam et al., 1986), since 
the precession motion of the rotor center with respect to the housing-centered position (zero static 
eccentricity) is assumed to be restricted to a circular orbit. With this approach the whole computa- 
tional domain is simply returning back to its quasi-steady-state condition, and the computational 
requirements are drastically simplified (see also fig. 28 for numerical gridding). This approach 
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Rj dty dz 


Figure 27. — Seal or bearing configuration in a 
constant axial plane Z. (Tam et al., 1988.) 



Bearing 
or seal 

z 


Figure 28. — Shaft housing configuration and fluid flow grid distribution. 
(Tam et al., 1988.) 


provides a very effective and competitive advantage enabling us to perform a true two-way, solid- 
fluid coupling analysis with the inclusion of a solid dynamic model simultaneously interacting with 
the fluid flow. Part 1 of appendix A highlights the governing equations formulated, the influence of 
turbulence, and the numerical algorithm employed on the multidimensional modeling of whirling 
seals and bearings. 
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Analytical Method of Bently and Muszynska (Lumped-Parameter Model) 

The lumped-parameter-model fluid average circumferential velocity ratio 2 X (Bently et al., 
1986; and Muszynska, 1986a,b) describes the fluid force more accurately than traditional analytical 
bearing and seal coefficient modeling. In this model the axial flow is uncoupled from the 
circumferential flow and affects the circumferential flow in a parametric way only. The model 
directly correlates to rotation-related circumferential flow parameters that are major factors in rotor 
stability. 

The lumped-parameter analytical modeling (see part 2 of appendix A) is divided into two parts. 
The first part deals with the fluid film forces within the seal or bearing. The second part couples the 
component with the rotor to account for two complex modes. The fluid film analysis demonstrates 
the importance of X to seal and bearing stability (Muszynska, 2001). The rotor/fluid-coupled mode 
solution is outlined in appendix A, and to a first order the coefficients of dynamic stiffness ( K = K D 
+ j Kq) become 


^ 0 = direct dynamic stiffness = Rotor DDS +Fluid DDS 




2 

Kq —M j (co p — Xco p ) 


fluid 


( 5 ) 


Kq = quadrature dynamic stiffness = RotorQ^s + Fluidqpjs 

= [^Mrotor +[ Z) ( ( °P - Xo) fl)]fiuid (6) 

Equations (5) and (6) illustrate that system stability can be readily altered (1) by controlling swirl 
(von Pragenau, 1982, 1985; Muszynska et al., 1988b; and Muszynska and Bently, 1989) and (2) by 
increasing radial fluid film stiffness (e.g., hydrostatic and ambient pressures; Bently et al., 2002). 
With control of both issues, injection and pressure, rotordynamic systems can be “walked through” 
rotor criticals (see Bently et al., 2002) during both turbomachine runup (powered) and rundown 
(unpowered) (also consider throttle chops). The rundown is most dangerous because the turbomachine 
is unloaded. 

The major effort is then to define the effects of X, fluid film stiffness, and damping. The rotor 
and fluid direct (real) K D and quadrature (imaginary) Kq components follow from the nonsynchronous 
forward perturbation numerical results (figs. 29 and 30), which show a large fluid inertia coefficient. 
When the fluid film component is coupled with the system and fluid inertia is neglected, stability 
thresholds agree with experimental evidence. 


Modeling Results 

Computational results. — Three-dimensional-flow numerical results show significant changes 
in the local values of the fluid dynamic forces along the seal (or bearing). They also show the 
existence of significant secondary flows and local separations even in cases of a large axial pressure 
drop. 


2 Fluid average circumferential velocity ratio multiplied by rotor speed X(ti R is the speed at which the fluid force rotates, 
where X can be identified as the ratio of the self-excited fluid whirl frequency to the rotational speed. By using hammer 
impulse testing of the rotor at constant speed, X(S) R can also be identified as a natural frequency of the rotor and fluid 
system. In many cases the value of X is near 1/2, but it can be greater or less than that depending on several factors (e.g., 
pressure gradient or viscous or injection effects). 
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3x106 



Figure 29. — Quadrature dynamic stiffness versus 
perturbation frequency. Identification of fluid 
radial damping D (8860 kg/s) as slope and velocity 
ratio \ (0.1 25) from the ordinate-intercept (-DX cop) 
of the straight line. Amplitude of rotor precession, 
A = c/2; change in fluid pressure, A p, 17 bars. 

(Tam et al., 1988.) 


Case 2 Rotational 



Modified perturbation frequency, 

(cop - \ cop) 2 , rpm 2 

Figure 30. — Direct dynamic stiffness as function of modified 
perturbation frequency. In these coordinates the results 
represent a straight line for easy identification of fluid 
radial stiffness K (intercept, 2.04 xIO 6 N/m) and fluid inertia 
M (slope, 4.63 kg). Fluid average circumferential velocity 
ratio, 0.125; amplitude of rotor precession, A = c/2; 
change in fluid pressure, A p, 17 bars. (Tam et al., 1988.) 


The integrally averaged forms of fluid direct K D and quadrature Kq dynamic stiffnesses arise 
from fluid radial F r and tangential F f force components divided by precession amplitude A. The 
dynamic stiffnesses are plotted versus perturbation frequency C0p (fig. 29) and in its modified form 
(fig. 30), where only forward perturbation was considered (C0p and CDp in the same direction). (See 
also part 1 of appendix A.) 
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Figure 31 . — Fluid (a) direct and (b) quadrature dynamic stiffnesses as function 
of rotational speed for two axial Reynolds numbers. Amplitude of rotor 
precession, A = 0.24c. A comparison of calculations and measurements 
by Childs in NASA CP-2238 and contract NAS8-33716. (Tam et al., 1988.) 


Comparison with experimental data. — Figures 31(a) and (b) compare numerical results 3 with 
the experimental results obtained by Childs (in NASA CP-2250, 1982, and in NASA CP-2338, 
1984). Although there were three steps within the length of the experimental seal, for simplicity the 
seal modeled was considered as a straight cylindrical seal without inlet preswirl. The experiment and 
the calculations qualitatively agreed. The calculated direct dynamic stiffness decreased with speed 
but did not exhibit the negative values of the experimental data at elevated rotor speed. The omission 
of the internal steps may account for this discrepancy. The quadrature dynamic stiffness increased 
with rotational speed and agreed reasonably well with the data of Childs. 

Nonlinear considerations. — For amplitude A = 0.5c (where c is bearing or seal radial clearance), 
the dynamic stiffnesses are higher than for A = 0.24c. The fluid radial stiffness K is a nonlinear 
increasing function of shaft eccentricity (the most pronounced effect), and pressure and fluid average 
circumferential velocity ratio X is a decreasing function of eccentricity (lesser effect). The effect of 
eccentricity on the fluid radial damping D and the fluid inertia in the range of parameters considered 
is not very high. 

Seal secondary flow. — Lower rotational speed and higher perturbation speed intensified the 
secondary flow zones and, by lowering the average circumferential velocity, enhanced rotor 
stability. The zone of secondary flow depends on the inlet boundary conditions and rotates with rotor 
precession frequency C0p. 


3 Modeling fluid, bromotrifluoromethane, or Halon 1301 (CBrF 3 ), used for fire suppression. 
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Figure 32. — Direct dynamic stiffness as function of rotational speed squared (synchronous perturbation). 
M( 1 - \) 2 = 0.0648 kg; fluid radial stiffness coeffecient, K, 1 .98x1 0 6 N/m; pressure-related circumfer- 
ential average velocity ratio, \ e ) = 0.1 1 1 (1 - X); amplitude of rotor precession, A = 0.24c. Note influ- 
ence of preswirl on average circumferential velocity ratio as well as increase of fluid radial stiffness 
with pressure differential. (Tam et al., 1988.) 


In the numerical model the preswirl was imposed at the seal inlet in either the forward or 
backward circumferential direction in the following form: 

V ps = RiCO R [l-r/R] l/3 (7) 

where Vp S is the preswirl circumferential velocity at radius r and Rf and R are shaft and seal radii, 
respectively. Figures 32 and 33 show calculations for seal flow with preswirl at two values of axial 
pressure drop. The preswirl clearly modified X, here shown for synchronous perturbation (C0p = COp). 
For higher axial pressure drops the circumferential flow is reduced as shown schematically in 
figure 34. 

Figure 35 illustrates the influence of inlet preswirl on the secondary (backward) circumferential 
flow field. It is evident that the preswirl intensified the recirculation zone if it was introduced in the 
direction opposite to rotation (lower resultant X and higher rotor stability). Preswirl in the direction 
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Figure 33. — Quadrature dynamic stiffness as function of rotational speed (synchronous perturbation). 
Fluid effective damping, D(1-\) = 8240 kg/s; pressure-related circumferential average velocity ratio, 

\ e , 0.141 (1-X); fluid tangential stiffness coefficient, K xy , 2.58 xIO 5 N/m; amplitude of rotor precession, 
A = 0.24c. Note influence of preswirl on average circumferential velocity ratio as well as increase of 
fluid radial stiffness with pressure differential. (Tam et al., 1988.) 



Figure 34. — Interpretation of axial and circumferential flow resulting in "spiral" flow, 
(a) High pressure drop, (b) Low pressure drop. (Tam et al., 1988.) 
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Axial flow 



Figure 35. — Influence of preswirl on circumferential flow 
field, (a) No preswirl, (b) Preswirl against direction of 
rotation, (c) Preswirl in direction of rotation. Shadowed 
areas indicate secondary flows. (Tam et al., 1988.) 


of rotation weakened the recirculation zone and intensified forward circumferential flow (higher 
resultant X). 

Effects of surface curvature and influence of rotating flows on turbulence. — As cited in part 1 
of appendix A, in the regions of negative velocity gradients the Richardson flux is also negative and 
a smaller value of 8 decay intensifies the turbulence level. Negative Richardson flux increases from 
the shaft walls (from inner to outer walls) and becomes a maximum at the separation regions. Then, 
the gradient of Richardson flux drops sharply and changes to a positive number where the turbulence 
level is damped and a secondary recirculating zone appears. 

In the case of forced swirl, formed by a spinning shaft, the centrifugal forces near the inner wall 
support the growth of vortical structures more than in curved-channel flow. The angular velocity 
gradients are generally negative all across the flow path, with the largest absolute values near the 
housing, where the largest turbulence will be generated. This situation is inherently unstable because 
high turbulence levels suppress the formation of secondary flow patterns. However, by introducing 
preswirl or multiple injection on the seal or bearing, often a counterswirling vortex can be formed 
near the housing wall, which is known in practice to have a stabilizing effect on rotordynamics. In 
this case the largest angular velocity gradient is shifted toward the shaft. Within the countervortex 
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Injection 

With rotation 7 = arctan (3) 

Radial 7 = 0 

Against rotation 7 = arctan (-3) 



Figure 36. — Injection geometry. (Tam et al., 1988.) 



Rotational speed, u > r , rad/s 



Rotational speed, wr , rad/s 


Figure 37. — Direct (a) and quadrature (b) dynamic stiffnesses as function of 
rotational speed for cases with fluid injection (synchronous perturbation). 
Change in fluid pressure, A p, 0.5 bar. (Tam et al., 1988.) 


zone both the Richardson flux and the gradient of angular velocity change sign, in turn reducing the 
turbulence levels in that region and enhancing system stability. 

Fluid injection effects . — Figures 36 and 37 illustrate the injection geometry and the influence 
of fluid injection on dynamic stiffness and give the corresponding angular relations. Three cases with 
injection were considered (radial or “vertical”, in the direction of rotation, and against rotation). 
Injection rate was proportional to rotational speed. All parameters exhibited nonlinear character as 


NASA/TM— 2004-21 1991/P ART2 


39 


functions of rotational speed. Fluid radial stiffness increased, fluid inertia decreased, and fluid 
effective damping D{ 1 - X + X t ) increased when injection was against rotation (X i < 0, where X i is the 
injection-related average circumferential velocity ratio) and decreased when injection was in the 
direction of rotation ( X • > 0). Radial injection caused a higher radial stiffness K and a lower X than 
in the no-injection case, a very important observation for pressurized, fluid-lubricated bearings as 
well as for antiswirl seal techniques. The average axial leakage was reduced by 7 percent. (See also 
part 2 of appendix A.) 

In developing the LE-7A liquid hydrogen-liquid oxygen engine for the H-IIA launch vehicle, 
subsynchronous whirl at 250 Hz (rotations/s) (15 000 rpm; 1571 rad/s) was found in the fuel 
turbopump. The instability could be controlled by tightening the bearing clearances, but that solution 
proved too costly. A CFD analysis was undertaken to assess the effect of Thomas forces (Alford 
forces) (fig. 38). The computational sector model included the power-stream flow path through the 
nozzle and one turbine stage of the high-pressure fuel turbopump (inlet pressure, 21 .4 MPa at 750 K) 
rotating at 703 Hz (rotations/s) (42 200 rpm; 4419 rad/s) and developing 28 800 kW of power at 
design conditions. Figures 39 and 40 give streamlines and pressure distributions at three normalized 
tip clearances (0.5,1 .0, 1.5). With these results rotor torque can be evaluated. Seal swirl or fluid-force 
circumferential velocity X causes the rotor system to become unstable. Reversing inlet swirl 
stabilizes the system as also noted by Thomas (1958), Benckert and Wachter (1978), and Tam et al. 
(1988). Altering the injection of liquid hydrogen coolant flows (seal leakage) at the seal inlet from 
normal to the rotor to tangential to the rotor mitigated the instability (figs. 41 and 42). 

In terms of Thomas forces as described by Motoi et al. (2003), the destabilizing force F f 
illustrated in figures 27 and 38 becomes 


s-X 


R 

2 R 


1 +- 


-T 0 


COS0 


0<6<2k 


a 


_7j-7b 


( 8 ) 


where l\ is the torque at maximum tip clearance, Y (1 the torque at average tip clearance, and '/’ uua | 
the total torque of all turbine blades. 


Kxy = F/e ( 9 ) 

where 8 is the ratio of eccentricity e to clearance c. As suggested by Thomas (1958) 

= (10) 

where L is the blade length, D the turbine diameter, and 

(3 = oc/(28 RF t ) (11) 

was evaluated though computation and agreed with experimental data. For the turbopump configu- 
ration of Motoi et al. (2003) K xy = 8.1 X 10 5 N/m and |3 = 0.56. 

The computational effort represents a more realistic approach by including nozzle and blade 
flows at three eccentric geometric conditions. For these rotor and stator computations, only steady 
solutions were considered, and “boundary data were averaged” and then passed off to the other 
subdomain at this boundary. The unsteady flow field can be computed but requires very large 
computational times, mesh sliding, and time resolution. Although the model is simplified with respect 
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Figure 38. — Thomas force. (Motoi et al., 2003.) Courtesy ISCORMA-2. 



(a) (b) (c) 

Figure 39. — Streamlines and pressure distributions, (a) Tip clearance, 0.5. (b) Tip clearance, 1 .0. 
(c) Tip clearance, 1.5. (Motoi et al., 2003.) Courtesy ISCORMA-2. 



Figure 40. — Effect of tip clearance on pressure distribution, (a) Tip clearance, 0.5. (b) Tip clearance, 
1.0. (c) Tip clearance, 1.5. (Motoi et al., 2003.) Courtesy ISCORMA-2. 
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Figure 41. — Damping ratio of each mode, (a) Turbine seal (cross section 
of fuel turbopump), (b) Conventional type, (c) Improvement type. 
(Motoi et al., 2003.) Courtesy ISCORMA-2. 



(a) (b) 

Figure 42. — Displacement in radial direction of fuel turbopump at engine firing tests, 
(a) Conventional seal, (b) Improved seal. (Motoi et al., 2003.) Courtesy ISCORMA-2. 
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to nozzle and blade time-dependent interactions and dynamic flow interactions between the casing 
and the rotor blade tip, Motoi et al. (2003) cite good agreement with data. Further, using these 
computational results as guidelines — instead of changing the turbine design to reduce the Thomas 
(Alford) forces and correcting the inlet swirl, which is known to stabilize seals (e.g., labyrinth seals) — 
provided the required turbopump stability. 

As another example of the effects of preswirl and injection, Kanki et al. (2003) investigated 
partial-admission turbine sealing. Gas turbine and steam turbine combined-cycle operation provides 
higher plant efficiencies at greatly reduced carbon dioxide, NO x , and other emissions than 
conventional systems. Partial-admission turbines are used to control load demand and therefore 
engender complex effects that often lead to subsynchronous vibrations. Under these conditions inlet 
swirl produces similar destabilizing effects in both cylindrical and labyrinth seals. Because the 
destabilizing force in partial-admission turbines depends strongly on swirl distribution, adding a 
swirl break (fig. 43) had a large effect on both full- and partial-admission turbines. Figure 44 gives 
natural frequency and damping ratio data for labyrinth and flat annular (smooth) seals and shows 
peak values for the smooth seal nearly twice those for the labyrinth seal with no swirl break inserts. 


M 




Figure 43. — Configuration of test seal, (a) Annular seal, (b) Labyrinth seal, (c) Flow 
blockers. Dimensions in millimeters. (Kanki et al., 2003.) Courtesy ISCORMA-2. 


NASA/TM— 2004-21 1991/P ART2 


43 


15 


N 

I 

S' 10 

0 

3 

cr 

0 


5 5 

3 

CO 


(a) 


A A A 

r ^ ♦* *: 



i : : • 

4/4 Admission . 

:• ; ; a 

3/4 

;* ; ; A 

2/4 

. .... X 

1/4 


12 3 4 

Seal pressure drop, A P, kPa 


0.05 


0.04 

0 

^ 0.03 

O) 

c 

1 0.02 

0 

Q 

0.01 


0 

Seal pressure drop, AP, kPa 


m *- ; 



u ;• ; 

• 

4/4 Admission 

pc 

. . . A 

3/4 



▲ 

2/4 

;• • • A- ;-A 

X 

1/4 

:• :-xa • 



: -a- • - 

• • A- 



A; • - 

; !■ 

: : x 

- -A - - A- 

■ A# 

.A 


; ; K 



< ^4-^ 



0 12 3 


N 

I 

0 

3 

O’ 

0 


2 

+- « 
0 


(b) 











4 * 






... 




1 

































• 

4/4 Admission . 

3/4 

2/4 

1/4 












A 




X 

i ! 




Seal pressure drop, AP, kPa 


0.05 


0.04 
^ 0.03 

O) 

c 

§- 0.02 

0 


0.01 


Seal pressure drop, AP, kPa 







i 4/4 Admission 

> 3/4 
^ 2/4 
< 1/4 





i 

... L 





i 





) 

< > XA 

II -Xi 






L*A- - -t 

■ l 





















A 

A 





> 

* X' ' 

< 

A A A g 

‘ A-A — * 





Figure 44. — Comparison of flat annular (smooth) seal and labyrinth seal, (a) Smooth seal, 
(b) Labyrinth seal. (Kanki et al., 2003.) Courtesy ISCORMA-2. 


For comparison, figure 45 shows cases with inlet swirl breakers only and with multiple-chamber 
swirl breakers (three labyrinth chambers and inlet). 

Figure 46(a) gives average swirl velocities and pressure drop data for four 90° flow blockage 
sector configurations (e.g., four of four sectors open (two transducers per sector; i.e., a full-admission 
condition)), and figure 46(b) gives data for three of four sectors blocked (i.e., a one-quarter- 
admission turbine). Data are also given for three of four and two of four sectors blocked. 

Conditions and parameters of Kanki et al. (2003 ). — Kanki et al. (2003) studied six-tooth 
labyrinth and smooth seals under the following conditions: experimental seal length (50 mm), 
diameter (129.5 mm), radial clearance (0.25 mm), casing surface (smooth), support shaft diameter 
(5 mm), support span length (300 mm), casing thickness (10 mm), inlet air pressure (0.4 MPa), outlet 
pressure (atmospheric), maximum flow rate (2.5 N-m 3 /min or standard m 3 /min), and rotational 
speed (zero). The pres wirier was a set of 45° helical gears with a standard 2-mm involute profile. The 
labyrinth tooth was 1 mm at the casing wall by 1 mm high and had a 0.5 -mm forward facing (surface) 
at the rotor interface. The cavity length was 8 mm with 2-mm inlet and exit lengths. The flow blocks 
were rubber inserts 1 by 3 by 8 mm for the labyrinth chamber or 1 by 3 by 2 mm for the inlet. Blocks 
were inserted around the circumference to block circumferential flow in either one, three, or five 
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Figure 45. — Effect of swirl breakers, (a) Inlet swirl breaker only, (b) Multiple-chamber swirl breakers 
(three chambers and inlet). (Kanki et al., 2003.) Courtesy ISCORMA-2. 


labyrinth chambers, with four blocks per chamber almost coinciding with the admission flow blocks. 
The damping ratio was derived from the free decay vibration caused by an impulse hammer striking 
the rotor. The rotor speed did not alter the results for this model and low-speed conditions. The 
computations were carried out using the theory of Kostyuk (1972) on a 130-mm-diameter, 0.25-mm- 
radial-clearance labyrinth seal with 6- and 9-mm cavity chamber height and width, respectively, and 
zero rotational speed. 

These examples illustrate the conflicting demands of negative preswirl for rotor stability and positive 
preswirl to mitigate turbine cavity heating and flow blockage losses of turbine blade cooling fluid. These 
conflicts must be properly balanced and resolved early in the design stage (see also part 3). 

Effects of axial flows . — Figure 47 illustrates a rotor housing study with large clearances to 
determine fluid effects in the convergent and divergent zones. At the smallest clearance the axial flow 
velocity was consistently zero, but the oil was axially pumped into and out of the divergent and 
convergent zones, respectively, with secondary flow zones in the high-clearance zone. In flow 
visualization experiments these zones were clearly observed (Braun et al., 1987). However, for this 
case these patterns were significantly changed by small axial pressure gradients (e.g., 0.1 bar; 
1.45 psi), and secondary flow zones tended to be small or vanished altogether. 
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Figure 46. — Measured swirl velocity, (a) 4/4 admission, (b) 1/4 admission, 
(Kanki et al., 2003.) Courtesy ISCORMA-2. 
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Case 7 


Axial flow velocity, V a , m/s 




Figure 47. — Mean axial flow distribution along unwrapped 
circumferential direction and secondary flow patterns 
of a bearing at various axial planes. Pressure drop, 0; 
rotor precession frequency, w p , 0. 


Summary 

Although advances in computational methods have increased our capability to determine forces 
in complex systems, feedback and validation of field data remain primarily in the proprietary realm, 
complicating information exchange. In one such exchange Thompson provided an average damping 
parameter for 1 15 individual compressor stages contained in 20 multistage bodies, but details are not 
likely to surface. 

Some physical descriptions of the instabilities given by, for example, Armstrong seem classic. 
Armstrong, in commenting on Alford’s paper, writes that at a certain temperature the frequency of 
mode vibration in an annular air cavity with six wavelengths is equal to the sixth diametrical mode 
of the stator. Under these conditions air gap variations in the labyrinth due to the stator vibration 
excite vibrations in the cavity. In turn, the resulting cavity pressure oscillations act on the stator and, 
if the phase is correct, self-excited stator oscillations occur. The mechanism could be stopped by 
inserting 24 equally spaced radial baffles across the air cavity. 

The stability work and rules of thumb of Alford and Thomas on turbomachine drum and shroud 
seals continue to be designed into aeronautical and aeronautical-derivative gas turbines and many 
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industrial steam turbines. Computational fluid dynamics (CFD) coupled with aeroacoustics and 
aeroelasticity are potential methods for validating these rules of thumb and have been developed. 
However, the field validation data and computational time associated with such a systematic program 
are unlikely to be supported unless a costly failure occurs. For example, the classic “fix” is still to 
add a split (severed) damper ring or cylindrical damper drum as a cure for aeronautical engine seal 
failures. 

Compressor blade tip vortices can engender a local stall region that grows into a rotating stall cell. 
Stage whirl instability is instigated when a blade stalls, developing a separation region, unloading 
the blade, and failing to produce sufficient pressure to maintain flow around it. The differential angles 
of attack cause the stall zones to propagate from the pressure side to the suction side of each blade 
opposite the direction of rotation. Storace et al. demonstrated the predominance of backward whirl 
over the operating range of a statically offset rotor in a four-stage, low-speed compressor. They also 
found that doubling the rotor tip clearance dropped torque efficiency by 1 .6 points versus 0.9 point 
for doubling the hub stator clearance. 

The data of Childs and field data reported in NASA Conference Proceedings served as sources 
of design information for the development of computational codes by Schuck and Nordmann, San 
Andres, and Tam to deal with instabilities. Equally important are the data from Bently Rotordynamics 
Research and Development Corporation and now Bently Pressurized Bearing Company. 

The rotordynamics codes of Nordmann, INDSEAL developed by Shapiro (1996), SCISEAL 
developed by Athavale, and TAMSEAL developed by San Andres, the three-dimensional fluid 
dynamics work of Tam et al., the CFD model, and the “lumped”-fluid-parameter model of 
Muszynska provide a foundation for studying the dynamics of rotor seal and bearing systems. (See 
also Hendricks, 1996.) 

These codes, the work of Tam et al., and the Bently and Muszynska models related herein 
corroborate and implicate that destabilizing factors are related though fluid average circumferential 
velocity. However, the dynamics of the rotational flow field are complex with bifurcated secondary 
flows that significantly alter the circumferential velocity. Several numerical experiments were 
carried out with the Tam et al. and CFD models and Bently and Muszynska’ s lumped-parameter 
model. Fluid injection and/or preswirl of the flow field opposing the shaft rotation significantly 
intensified these secondary recirculation zones and thus reduced the average circumferential 
velocity, whereas injection or preswirl in the direction of rotation significantly weakened these 
zones. A decrease in average circumferential velocity was related to an increase in the strength of the 
recirculation zones and thereby promoted stability. The lumped-parameter model describes the 
dynamic forces for relatively large but limited ranges of parameters and is extremely useful in system 
dynamic analysis of rotor bearings and seals. 

Seals play a prominent role in stabilizing turbomachines, and understanding of sealing forces in 
high-power-density machines has emerged, although the details continue to be analyzed yet today. 
Correcting the seals enabled the space shuttle’s high-pressure fuel turbopump to operate success- 
fully. Also, correcting the fluid injection to control seal preswirl of the fuel turbopump enabled the 
LH-7A hydrogen-oxygen engine to operate safely. 

Finally, if the engine and power conversion industry wants a smooth-operating, reliable, long- 
lived machine, it must pay very close attention to sealing dynamics. Does this mean that gas path or 
power-stream dynamics are not important? Of course not. Brayton cycle machines, for example, are 
fairly mature. However, turbomachine sealing dynamics represents your competitive edge. 
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Appendix A 

Numerical and Analytical Modeling 

Part 1 — Numerical Modeling 

For an observer located in the rotating frame the relation between absolute, relative, and grid 
velocity becomes 


Ugfos—U+QdpXr (A— 1) 

where j/ abs is the fluid absolute velocity, fj the fluid relative velocity, c5p the rotor precession 

(perturbation) frequency, and r the local radius in the fluid. For any given perturbation speed ratio 
the transformed moving-coordinate grid system must be consistent with the specified rotor preces- 
sional speed imposed (fig. 27). 

Governing equations . — Conservation of mass can be written as 

|-(p/)+V-(pt/V=0 (A- 2) 

where p is the fluid density and J the determinant of the Jacobian matrix (i.e., the physical volume 
of a grid cell in the finite difference approach). 

From the Euler identity and the small circular orbit assumption the rate change of the Jacobian 
matrix can be made to vanish. The “quasi-steady” analysis is valid and equation 
(A- 2) takes the final form 



p u u=o 


(A-3) 


The three-dimensional, turbulent, time-averaged, “conservative” 4 Navier-Stokes (N-S) equa- 
tion in the coordinate system rotating at the perturbation angular velocity (0 p is 


5( p UJ 

~dT 


+V- p JUU \=-JVp-J 


Centrifugal 

' I ’ 

ptDp x(c5p xr )+2p(0p xU 

T 


Shear 

+/Vx+Jpg 

T 


04-4) 


Coriolis Body force 


It is very important to note the centrifugal and Coriolis source terms arising from the shaft whirl- 
ing at (Op. 


Conservation is respected at every grid point in the mesh; viscous fluctuating terms are time averaged. 
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In this equation the shear stress tensor Vx contains time-averaged Reynolds stress terms. The 

turbulence model used in our seal flow simulations makes use of the eddy viscosity concept to 
compute those terms, resulting in the following expressions: 


7Vx = -pUjiijJ = \i t 


dui + ^ u j 
dxj dxf 


-V-t/8« 
3 J 


(A- 5a) 


where the turbulent eddy viscosity is proportional to the mixture density of the turbulent flow p, 
the rate of strain (]), and the square of the mixing length of turbulence l m . The centrifugal and Coriolis 
forces appear to be caused by the coordinate rotation and are incorporated into the computational 
model as additional volumetric momentum sources. 


S(u) = -2pC0pV Vol ce n, S(v) = |pc0pr + 2pC0pwj Vol ce u (A -5b) 

where Vol ce p is the physical volume of a grid cell (from the Jacobian J). 

Turbulence modeling for seal and bearing flows . — The characteristics of seal and bearing 
flows can be summarized as strong swirl, high viscosity, and high pressure drop associated with 
significant wall roughness and small clearances. Several attempts at predicting rotating turbulent 
flows indicate that using complex full Reynolds stress equations was only as good as using those 
obtained with simple correction of the turbulence generation and dissipation k- 8 . Furthermore, 
existing turbulence models can successfully predict only simple boundary layer or mildly recircu- 
lating flows. 

Effects of surface curvature and rotating flows on turbulence : Several investigators have 
attempted to extend both the algebraic Reynolds stress and k-e models for turbulent rotating flows. 
Salient features of the effects of surface curvature and rotating flows on turbulence modeling can be 
discussed briefly. Most relevant approaches are based on Bradshaw’s suggestion (Bradshaw, 1973. 
See also Hoffman etal., 1985; Muck etal., 1985;andSmits etal., 1979.) that the streamline curvature 
effects can be modeled by modifying the turbulence length scale. The two most relevant approaches 
are those of Rodi (1979) and Launder et al. (1977). Both incorporate the “rotating” correction into 
the constants of the turbulence dissipation rate term. The extension of Rodi introduces a correction 
to the model constant in the generation part of the 8 equation source term. The Launder extension 
modifies the destruction part of that term in a similar manner by introducing the Richardson flux 
number: 


D _ k 2 u dru 
£ r or 

Both the Rodi and Launder et al. approaches respond to the angular speed, to the swirl velocity 
gradient, and to the turbulence time scale. Analysis of the Richardson flux correlation indicates that 
in the regions of negative velocity gradients R is also negative and smaller 8 decay intensifies the 
turbulence level. Negative Richardson flux starts increasing from the shaft walls (from inner to outer 
walls) and becomes maximum at the separation regions. Then, the gradient of Richardson flux drops 
sharply and changes to a positive number where the turbulence level is damped and a secondary 
recirculating zone appears. 

Two-equation model extensions: In extending the two-equation turbulence model, like other 
turbulence flow researchers (Rodi, 1979, and Bradshaw, 1969, in particular), Kim and Rhode (1999) 
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have also incorporated Richardson flux corrections, modified the turbulence mixing-length scale, 
and elaborated on the log-law wall function with an attempt to account for swirling streamline- 
curvature flow influences on the k-e turbulence model. They claimed only very limited success. 

For coupled swirling flows these models fail to capture the essential physics of usually 
anisotropic, nonhomogeneous, and three-dimensional turbulence structures. However, the general 
consensus is that the empirical constants in the semi-empirical k- 8 model are invalid for complex 
flows anyway. For these reasons a simple Prandtl mixing-length model was used to represent the 
turbulence in seal and bearing flow passages. Further, the purpose herein is to develop the dynamics 
of fluid- solid-coupled modeling rather than to develop a turbulence model for nonplanar interfaces. 

Prandtl mixing-length turbulence model. The zero-equation Prandtl mixing-length model is 
simple and well established, and the mixing-length scale is determined by the clearance. The 
turbulence model used the eddy-viscosity concept 


where £ is the normal distance measured between the stator and rotor surfaces, y is the minimum local 
distance measured from the stator and rotor surfaces, and b and a are coefficients. 

Numerical algorithm , boundary conditions , and grid . — An implicit pres sure-correction method, 
the SIMPLEST algorithm proposed by Spalding (1981), was used to solve the numerical model. The 
momentum equations were first solved with a guessed pressure field. An implicit successive iterative 
method with a pressure-correction equation derived from the continuity equation was used in the 
algorithm. Subsequent pressure and momentum were corrected to satisfy continuity. The whole 
process was repeated until satisfactory convergence was achieved. 

At the flow inlet both the pressure and the circumferential velocity (with or without preswirl) 
were specified. Only pressure was specified at the outlet. For a specific perturbation speed ratio the 
proper rotor rotational speed and the rotor precessional speed were specified. The transformed 
moving coordinate grid system must be consistent with the specified rotor precessional speed 
imposed. The no-slip velocity condition at the rotor and stator walls and fluid shear stresses acting 
on the rotor and stator surfaces were generally incorporated. Shear stresses on the rotor and stator 
surfaces were calculated on the basis of the Couette flow assumption with universal “log-law” of the 
wall function approach adopted. For the roughness of the wall boundary specifications Childs et al. 
(in NASA CP-2133, 1980, and in NASA CP-2443, 1986) demonstrated that surface roughness 
depends on the flow direction (i.e., there is higher resistance to flow in the axial than in the 
circumferential direction). Consequently, different rough-to- smooth correlations were specified in 
the axial and circumferential directions (the coefficients 2.9 and 1.1 were multipliers of the smooth 
friction relations). Finally, cyclic boundary conditions in the circumferential direction were adopted. 

A three-dimensional, nonorthogonal, body-fitted computational grid has been selected for the 
computations (fig. 28). The grid allows for steep variations in physical parameters (fluid preswirl in 
particular). It has been also optimized according to a surface Reynolds number criterion (Re + < 100) 
of the imposed wall function to ensure proper resolution of the boundary layer near the rotor and 


Dynamic forces. — Fluid circumferential force components F r (radial) and F t (tangential) were 
calculated from the pressure distribution 


where R[ is the shaft radius, / denotes bearing or seal length, and 0 is the angle measured (clockwise) 
from the position of minimum clearance. 



a = 0.09 


stator walls. 



(A-6) 
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The third component of the fluid force F^is related to the axial flow such as that occurring in seals 
of pumps and compressors due to balance piston loading. In bearings F is usually small relative to 
F r and F p but for seals in machines of the SSME category the axial flow and pressure drop are large 
and F z has a significant value. However, it was assumed to remain implicitly linked to F f and is a 
source for future work. 


Part 2 — Analytical Lumped-Parameter Model 

The lumped-parameter model represents an analytical approach to rotordynamics with fluid film 
components (Bendy etal., 1986; and Muszynska, 1986, 1988a, 2001). The fluid force model was identified 
by using extensive modal perturbation testing (Muszynska and Bendy, 1990; Muszynska, 1995). 

The most important novel feature in the model is the fundamental assumption that the fluid 
dynamic force is related to the circumferential flow, which is generated by rotor rotation in the rotor/ 
stationary clearance. Note that this model is applicable for fluid-lubricated bearings, seals, and any 
other radial or axial clearances between a rotating and stationary part. 

At the beginning it is assumed that the rotor is centered within the clearance and the averaged 
(lumped) fluid force is rotating at the angular velocity X(ti R (not a constant co^/2, as is often assumed, 
e.g., Black and Jenssen, 1969-70), where co^ is rotor rotational speed. The fluid force expressed in 
coordinates rotating at angular velocity Xod R has three classical components: stiffness K (stiffness 
force proportional to shaft radial displacement z r ), damping D (damping force proportional to shaft 
radial velocity dzjdt ), and fluid inertia M (fluid inertia force proportional to radial acceleration 
d 2 zjdt 2 , where z r = x y + jy r is the complex displacement of the rotor within the clearance. In 
stationary coordinates the fluid force becomes 


- F = Kz + D(z - jX(0 R z) + m[z - 2jX($ R z - X 2 co \ z^ 
z = X + jy, j = V-L, \z\ = t]x 2 + y 2 , • = d / dt 


(A -7) 


where X is the fluid circumferential average velocity ratio (a measure of the circumferential flow 
strength), z is shaft radial complex displacement (conventionally, x denotes horizontal, andy denotes 
vertical), and z-z r exp (j) i(0 R ) is the transformation from stationary to rotating coordinates. 

In the bearing and seal coefficient format and in the chosen reference system (fig. 27) the fluid 
force for clockwise rotor rotation is as follows: 


1 
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~F t 
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(A -8) 
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Note that the cross-coupled stiffness term is generated here by the radial damping due to rotation and 
is proportional to X and CD^. The radial stiffness is modified by centrifugal inertia and carries the 
negative sign. “Cross damping” is due to the Coriolis inertia force. 
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When in response to external excitation by nonsynchronously rotating force P exp(/cOp0> the 
rotor precessional (orbital) motion is circular with amplitude A, phase a, and frequency C0p that 

is, 


z=Ae j((o P t+a) (A- 9 ) 

When equation (A- 9) is substituted into equation ( A-7 ), the following relation results: 


Pe j(0pt = Ae j((0pt+a) 


2 

K + /D(co p — X-CO p ) — M(co p — Xco ^ ) 


(A -10a) 


The external force divided by the rotor response yields the fluid-related complex dynamic stiffness 

=K-M{(0 P -X(0 R ) 2 +jD{(0 P -X(0 R )^K D +jK Q (A-lOb) 

Muszynska and Bendy (1990) and Muszynska (1995) identified the direct (real) K D and the quadrature 
(imaginary) Kq components of the fluid dynamic stiffness in equations (A-10) by using nonsynchronous 
forward and backward perturbation testing. Tam et al. (1988) also obtained the elements of the fluid 
dynamic stiffness from numerical calculations. (Selected results are illustrated in figures 29 to 33.) 
Note the high value of the fluid inertia effect, M- 4.63 kg, exceeds by about 400 times the fluid mass 
in the seal (fig. 30). 

For an isotropic shaft carrying a disk of mass M x and a journal of mass M 2 and supported by a 
fluid film bearing on one end and a pivoting rolling-element bearing on the other end, with shaft 
partial stiffness K x (left side) and K 2 (right side) (fig. 48), the governing equations become 

M{z\ + D s z\ + ( K x + K 2 )z x ~ K 2 Z 2 = ^ 2 r(0 i ex P 0%0 (A - 1 1) 


M 2 z 2 +Mf(z2 -^\z2) + \D + yf2\ 

X ( Z 2 ~ AcO ^ Z 2 ) + (^0 + ¥ l)^2 + ^3 z 2 +^ 2^2 “^ l)]-0 ( A - 12 ) 

Zi(t)=Xi(t)+j yi(t\ 1=1,2 (A-13) 

where D is fluid damping and D s is external damping; Kq and M^are fluid film stiffness and fluid 
inertia; m 2 and r are the mass and radius of unbalance, respectively; and \|/ 1 and \| f 2 are nonlinear 
stiffness and damping functions of the journal radial displacement I z 2 l, respectively. 

Muszynska (1985, 2001) gives and discusses the solution for the stability threshold of pure shaft 
rotation: 


(o| T =(l/X) 


\ M U 


+ 


K 2 (Kq+K 3-M 2 K 1 IM{) 
M x (K 2 +K 0 +K 3 -M 2 K x !M x ) 


(A-14) 


For co^ < co| T the rotor is stable; for co^ > co| T the rotor becomes unstable. Because in equation 

(A-14) the fluid inertia and the journal mass are small relative to the disk mass, and shaft stiffness 
K x dominates other stiffnesses being in sequence, the threshold of stability can be approximated as 
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Zi (t) Z 2 (t) Z-\ (t) z 2 (f) 




P e/“ f 



Figure 48. — Rotor and fluid film modeling, (a) Rotor and seal model, (b) Rotor and bearing model, 
(c) Rotor system accounting for two-complex-mode model. Nonsynchronous perturbation force 
P applied at the journal (or at the disk) as in equations (A-1) and (A-12). (Muszynska, 2001 .) 


co| T =(1/X) a /(^i/M 1 ) (A- 15) 

Note that the rotor instability thresholds (eqs. (A-14) and (A-15)) are inversely proportional to X. 

The rotor and fluid model (eqs. (A— 11) to (A-13)) also provides stability criteria for the 
unbalance-related synchronous response of the rotor (Muszynska, 1986). It also explicitly provides 
the postinstability-threshold rotor self-excited vibrations, fluid whirl and fluid whip, as two varieties 
of the same phenomenon (fig. 16). Note that in equations (A-1 1) and (A-12) it has been assumed that 
the journal was centered within the bearing clearance. Muszynska (1986b) solved the more realistic 
cases considering eccentric positions for the journal. Below there is an example on how the 
nonconcentric cases can be handled. 

Implementation of the fluid-force model (eq. (A— 8)) into another one-disk model of an isotropic 
rotor, operating in a fluid environment is as follows: 


M^+D.i+^z+Dli-ifX-X^IzDco^JJ+f/r+xi/idzI^PexpO-S) (A- 16) 

where M d and K d are rotor modal mass and stiffness, respectively, for the first lateral mode and P 
and 8 are the magnitude and angular orientation, respectively, of a constant radial force applied to 
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the rotor. At \z\ = 0 (concentric case) X n - X; then, with increasing eccentricity, X n decreases to zero, 
which means that the circumferential flow is no longer dominant. In this model rotor unbalance, fluid 
nonlinear damping, and fluid inertia are not included for clarity. The model contains, however, a 
nonlinear function of fluid circumferential average velocity ratio X n , which is a decreasing function 
of the rotor radial displacement Izl . This model reflects the fact that when the rotor is concentric within 
the clearance, the circumferential flow is the strongest. A rotor at higher eccentricity suppresses the 
circumferential flow (some secondary flows occur and most often the axial flow then becomes 
dominant). The model (eq. (A- 16)) shows the nonlinear effect of the fluid force when constant 
external force displaces the rotor from the concentric position in the fluid-filled clearance. 

The instability threshold for the linear part of the rotor model (eq. (A- 16)) is as follows: 


As previously (eqs. (A-14) and (A-15)) the instability threshold is inversely proportional to the fluid 
circumferential average velocity ratio X. Changing this ratio would directly affect the rotor instability 
threshold. One way to achieve this change is to move the rotor to higher eccentricity by applying 
radial force. There is also another way. It has been shown that external injection of the fluid into the 
seal clearance in the tangential direction, either the same or opposite to rotor rotation, changes X as 
well as the fluid radial stiffness: 


The new X depends on the direction of injection, so it can increase or decrease. The fluid radial 
stiffness always increases with any direction of injection because it creates additional pressure in the 
clearance. As can be seen in equation (A-17), if the fluid is injected in the direction opposite to rotor 
rotation, the rotor instability threshold may be significantly increased (Muszynska et al., 1988b; 
Muszynska and Bendy, 1989; and Muszynska, 2001): 


Equation (A-19) clearly shows that external tangential injection (1) alters rotor stability and (2) 
increases radial fluid film stiffness. With control of both issues, tangential injection and pressure, 
rotordynamic systems can be “walked through” rotor instabilities (Bendy et al., 2002). 

The major effort for the CFD analysis is then to define the effects of X, which is a measure of 
circumferential flow strength, fluid film stiffness, and damping with coupling to the shaft dynamics, 
and thus solve the problem of fluid-solid coupling interaction. 
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X+X(|z|)^ X+X(|z|)±X ext , K+V! -»A'+y 1 +£ ext 


(4-18) 



(A-19) 



Appendix B 


Turbomachine Fluid/Structural Interactions 

All turbomachines respond to unsteady thermal fluid dynamic processes, but the structural codes 
discussed here closely replicate only structural dynamic imbalance. The catch is that these dynamic 
imbalances are of mechanical origin (i.e., added torque through weight placement) and conform to 
structural responses. On the other hand, thermal fluid dynamic perturbations are driven by the 
inherent unsteadiness of the flow coupled through fluid/structural interactions to either reinforce or 
damp perturbations. Unless these factors are built into the engine codes, resonant or unstable engine 
phenomena will always appear that may necessitate engine rework or at best altering component 
design. Although expensive in terms of development time and testing cost, experimental validation 
of computations is a necessary part of certification. 


Full-Scale Systems 

Industrial gas turbine . — As turbomachines become more sophisticated, the need for full 
component and engine assembly test and analysis becomes acute. Although these analyses are 
expensive to set up and computationally intensive to carry out, whole-engine gas path modeling is 
being developed by several researchers (e.g., Numerical Propulsion System Simulation (NPSS) by 
NASA GRC program manager Joseph Veres), and structural analysis codes exist (e.g., MSC- 
NASTRAN). 

Surial and Kaushal (2003) applied MSC-NASTRAN superelements to a Rolls-Royce 
aeroderivative, three-spool gas turbine with (1) representing the engine casing, (2) the low-pressure 
turbine (LPT), (3) the intermediate-pressure turbine rotor (IPT), and (4) the high-pressure turbine 
(HPT). The casing superelement (0) was modeled by using quadrilateral (QUAD4) and three-node 
(TRIA3) linear shell elements (first order) and standard beam elements. The casing superelement 
model had eight modules: 

1. Front end and low-pressure case 

2. Intercompressor duct 

3. Plenum 

4. Aeronautical intercase 

5. Combustor 

6. High-pressure case 

7. IPT case 

8. LPT case 

Other superelements included (10), the LP rotor; (20), the IP rotor; (30), the HP rotor; (40), the engine 
coupling; (50), the pedestal shaft; (60), the generator coupling; and (70), the generator rotor. 
Figure 49 illustrates the computational mesh of the aeroderivative gas turbine with figure 50 showing 
the LP rotor, IP rotor, and HP rotor or spools. The drive trains were also modeled, and Guy an 
reduction was mapped each onto a common centerline model. The mode shapes were animated by 
using PATRAN to determine whirl directions. Modal frequency analysis was performed by using 
NASTRAN (SOL 1 1 1) with gyroscopic effects considered by using NASTRAN sssalterridgyroa.705. 
Pedestal damping included CDAMP elements. An unbalance load of 7.22X 10 -3 kg-m (10 oz.-in.) 
was applied at seven engine locations giving seven excitation and nine computed responses. The 
results were in good agreement for those points measured (table III). 
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Rear mount 
(both sides) 


Figure 49. — Finite element model and engine-to-ground mounts. (Surial and Kaushal, 
2003.) Courtesy ISCORMA-2. 




(b) 



G 




Figure 50. — Illustration of some Guyan points, (a) On LP rotor, (b) On IP rotor, 
(c) On HP rotor. (Surial and Kaushal, 2003.) Courtesy ISCORMA-2. 
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TABLE III.— RESULTS OF CRITICAL SPEED ANALYSIS 
AND MEASUREMENTS 


Mode 

number 

Frequency, 

Hz 

Description 


Experimental 

Computational 


92 


9.42 

Drive train torsion 

99 


21.55 

Drive train torsion 

104 


25.13 

Engine and generator (axial) coupling 

109 

33.6 

33.73 

LP rotor mixed whirl 

111 


37.04 

LP rotor and generator rotor axial 

115 

42.72 

41.67 

LP rotor mixed whirl 

117 

47.3 

44.58 

LP rotor forward whirl 

121 


50 56 

LP rotor reverse whirl 

122 

53.4 

51.28 

LP rotor pedestal shaft first conical mode 

123 

53.4 

52.58 

Pedestal shaft second conical mode 


TABLE IV, — NONDIMENSION AL NATURAL FREQUENCIES FOR THREE BLADE MODELS 


Blade model 

Nondimensional natural frequency 

Cantilevered blade without root cross section 

0.5062 

0.9002 

1.8062 

2.3128 

3.7909 

Cantilevered blade with root cross section 

0.4681 

0.8388 

1.7777 

2.0235 

3.0766 

Blade with disk sector 

0.3190 

0.4747 

0.9715 

1.2170 

1.7808 


Industrial gas-turbine-driven compressor . — Camatti et al. (2003a) relate instability problems 
in full-scale testing of gas-turbine-driven, back-to-back compressors with greater details, and 
Camatti et al. (2003b) give a second system example. Table IV and figure 5 1 give seal locations and 
types. The impeller slant-toothed labyrinth seal and the interstage honeycomb balance-piston seal 
were unstable. The final balance-piston, abradable shroud seal, teeth on rotor (fig. 51(c)) appeared 
to be stable. The impeller shroud labyrinth seals (fig. 51(a)) were replaced by labyrinth seals with 
staggered, multiple, axial slots machined into the intertooth cavities and referred to as web seals 
(fig. 52). These slots were effective in reducing the impeller shroud circumferential fluid force 
velocity Q!k. The honeycomb interstage seal (fig. 51(b)) was modified by adding convergent taper. 
The operating taper was a compromise between the need to increase frequency and the need for 
sufficient damping and was set at 0.06 mm (0.0024 in.) radial. However, because of pressure and 
thermal effects the machined taper was twice that, or 0.12 mm (0.0048 in.) radial; these changes 
nearly doubled the damping. Nevertheless, even with changes in the impeller seal (or web seal; 
fig. 51(a)) and seal tapering (honeycomb), operation with plugged radial-diffused fluid injection 
holes, or “shunts,” provided only marginal improvements. However, opening the shunts to the 
honeycomb seal inlet (fig. 51(d)) decreased gas swirl (i.e., reduced QX; see also figs. 36 and 37), and 
the system operated without a trace of instability, even to theoretical surge at 120 bar (120 percent 
operational). Opening the shunts perturbed diffuser and back-side impeller fluid dynamics and 
balance and decreased or nearly eliminated interstage- seal (fig. 51(b)) inlet gas swirl (i.e., reduced 
or eliminated QX; see also figs. 36 and 37). Recalling that alterations of the CDP seal changed flows 
throughout the engine (e.g., fig. 12 of part 1), shunt fluid injection may also alter impeller flow and 
pressure fields and thus enhance power-stream stability. (See also stability criteria in Fulton, 2003). 
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Figure 51. — Compressor cross-sectional drawing showing detail of rotor plus seals, (a) Impeller 
shroud labyrinth seal, (b) Honeycomb interstage seal, (c) Abradable seal, (d) Honeycomb 
interstage seal. (Camatti et al., 2003.) Courtesy ISCORMA-2. 




Figure 52. — Web seal, (a) Sketch of conceptual web seal, (b) Photograph of web seal or slotted 
labyrinth seal. (Camatti et al., 2003.) Courtesy ISCORMA-2. 


Figure 53 shows the nature of this labyrinth and honeycomb dual-seal concept. The two labyrinth 
teeth are in the extreme foreground, and between the labyrinth teeth and the honeycomb is a recessed 
annulus. Within the annulus are shunt holes that connect the high-pressure diffuser to the annulus. 
These holes can be plugged or open. When they are open, high-pressure diffuser gas flows radially 
through them into the annulus. Because this annulus contains the highest pressure in the seal, some 
gas flows toward the foreground, across the two labyrinth teeth, and beyond, flowing radially 
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Figure 53. — Honeycomb seal rub region with provisions for shunt. (Camatti et al., 2003.) 
Courtesy ISCORMA-2. 



Figure 54. — Aspirating seal GE90 test, (a) At shutdown phase, (b) At steady-state operation. 
(Tseng et al., 2002.) 


outward along the back of the impeller disk. However, most of the injected flow goes through the 
clearance between the smooth shaft and the honeycomb cells. 

In many respects this dual-seal configuration functions like the aspirating seal (fig. 54). In the 
aspirating seal the labyrinth tooth controls leakage flow at elevated clearances (fig. 54(a)) and forms 
a high-pressure cavity for the diffused-injected flow in the engine cavity at operating conditions 
(fig. 54(b)), where flow is controlled by the narrow gap in the face seal (Tseng et al., 2002). 
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Padded-finger seal fluid. — Altering the geometry of the basic finger seal concept, Braun et al. 
(2003) developed a new seal that combines the features of a self-acting shaft seal lift pad as an 
extension of the downstream fingers with an overlapping row of noncontact upstream fingers. The 
upstream and downstream fingers are further reinforced each by a single row of overlapping fingers 
(figs. 55 and 56). These overlapping fingers reduce the axial and radial flows along the compliant 
fingers that allow for radial motion of the shaft/seal interface. The lift pad rubs the surface until 
sufficient force is generated to lift the pad into a noncontact position. Leakage will occur between 
the noncontact overlapping fingers and the lift-pad extended fingers. These seals respond to both 
radial and axial shaft perturbations and some degree of misalignment with minimal hysteresis; the 
latter is a problem in conventional finger seals and some types of early brush seal designs. The 
importance of including coupled thermal/fluid/structural interactions makes the difference between 
a seal that works and one that does not. With this in mind, Braun et al. (2003) were able to determine 
flow patterns and provide seal dynamic response for the padded- finger seal by relating the cantilever 
or finger stiffness to fluid stiffness. 

Fluid film journal bearing. — Temis and Temis (2003) also cite the importance of fluid/ 
structural interactions in predicting load capacity and stability for conventional fluid film journal 
bearings. In their analysis the energy equation and conjugate heat transfer are neglected as are 
compressibility effects, which need further consideration. Nevertheless, consider a bronze steel- 
shaft journal bearing. Young’ s moduli are £ bronze = 0.75 X 10 5 MPa and £ steel = 1 . 15 X 10 5 MPa with 
a Poisson ratio of 0.3 for both. Shaft diameter is 30 mm, clearance ( c = 8) is 0.05 mm, and the ratio 
of turbine diameter to blade length DIL = 1. The shaft rotates at 30 000 rpm, with oil viscosity of 
0.02 X 10 -6 Pa-s at an eccentricity ratio £ = % = etc = 0.8. 

Figure 57 illustrates the loading, figure 58 the deformed fluid film journal bearing surfaces, and 
figure 59 the actual displacements. Note that these displacements have to be reflected to circumfer- 
ential angle cp = 0. With a compliant interface the load exerted by the fluid film on the bearing support 
tends to a finite value rather than an infinite value as eccentricity £ = %— >1 (fig. 60(a)) and the attitude 
angle 0 tends to be nonzero (fig. 60(b)). These compliant surfaces greatly affect stability margins, 
tending to lower stability with increasing bearing D/L (fig. 61). The singularity (e = % — >1) was not 
examined because it represents boundary-lubricated or dry-contact rub with significant heat release, 
another feature to be added to bearing analysis research and development. (See also the section Life 
and Reliability Issues in part 3.) 

Rotor impact loading. — In Navy aircraft slammed onto the deck by the arresting hook cable and 
other hard landings, the aircraft gas turbine’s rotor bearing undergoes abrupt changes in loading, 
termed “impact loading,” a common occurrence. Aircraft tires are seldom spun to speed prior to the 
impact loadings of landing and associated spin-up to speed. Runway tire marks and maintenance 
records attest to the brutal nature of hard impacts absorbed by landing gear struts and tire spin-up. 

The analysis is beyond our scope, but the similarities are not as we consider yet another class of 
fluid/structural interactions. Muszynska (2003) shows the effects of horizontal rotor drop onto 
retainer or load-sharing bearings due to the loss of magnetic bearing support to be skid and bounce. 
Skid represents a form of dry or boundary-lubricated contact for very short periods of time. The time 
rate of change of the angular coordinate position is modified by subtracting the term KCl (tangential 
radial stiffness coefficient times rotational speed). The K term is related to surface shear force or pA 
(kinetic friction times rotor radial load and surface speed differential). The energy release, surface 
deformation, and lubrication problems have yet to be considered, as is also the case of A = Mg (rotor 
mass times gravitational acceleration) rotor loading where other perturbations can significantly alter 
these results (e.g., a vertical rotor). Depending on the value of K the skid or bouncing stage leads the 
rotor either to stabilization or to a very dangerous mode of self-excited vibrations called dry whip. 
Muszynska (2003) gives requirements for retainer (backup) and load-sharing bearings. They must 
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Figure 55. — Typical dimensions (in inches) and 
appearance of padded finger. (Braun et al., 2003.) 


/"Stationary pare/ 


c 


Low-pressure 
cavity (zone 3); 
second row 
finger pad 


Connecting 


/--Padless 

finger 

High-pressure 
cavity (zone 1) 

Leakage path 
^ ^ -under pad 
(zone 2) 


Surface rotating into 
plane of paper 


Figure 56. — Boundary conditions for finger seal. 
Schematic cross section with two rows of 
padded low-pressure and padless high- 
pressure fingers and definition of pressure 
zones. (Braun et al., 2003.) 



Figure 57. — Section of fluid film journal bearing. 
(Temis and Temis, 2003.) Courtesy ISCORMA-2. 



Figure 58. — Fluid film bearing with deformed 
surfaces (1 -bearing surface; 2-shaft journal 
surface). (Temis and Temis, 2003.) 

Courtesy ISCORMA-2. 


provide machine support at rest and at all transient states between rests and provide stable operation 
as well as support in the event of magnetic bearing failure. In addition, these bearings must ensure 
safe, efficient, cost-effective machine operation. 

Although first-stage disk failure of the high-pressure turbine is the critical component in 
turbomachine failure, turbine and compressor blade fatigue failures also represent serious design 
challenges. Rzadkowski et al. (2003) investigated the natural frequencies of shrouded blade/disk/ 
shaft models of one, two, and three turbine stages. The turbine disk modeled had a 0.685-m outside 
diameter with a 0. 196-m bore and 144 tuned blades 0. 159 m long (fig. 62). Details and specification 
of the model, geometry, element interface connections, and normalization parameters were not 
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Carrying force, 


0.005 



Figure 59. — Bearing and shaft journal radial displace- 
ments as function of circumferential angle; §1 = bearing 
displacement, I 82 I = magnitude of shaft displacement, 
8 b = total displacement gap. (Temis and Temis, 2003.) 
Courtesy ISCORMA-2. 



Figure 60. — Bearing resultant carrying force (a) and attitude angle (b) as function of ratio of 
eccentricity to clearance. (Temis and Temis, 2003.) Courtesy ISCORMA-2. 



Ratio of turbine diameter to blade length, D/L 

Figure 61 . — Boundary between stable and unstable 
orbits for shaft journal motion in bearing. (Temis and 
Temis, 2003.) Courtesy ISCORMA-2. 
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Figure 62. — Dimensionless relative displacements of three shrouded 
blade disks placed on shaft illustrating modal shapes, (a) Inner 
disk motion, (b) Outer disk motion, (c) Nearly equivalent motion. 
(Rzadkowski et al., 2003.) 
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disclosed, making it impossible to replicate or apply the authors’ results directly without their 
consent. Note that detuning the blades enhances stability. The modes of the bladed disks are 
classified by using an analogy with axisymmetric modes, which are mainly characterized by nodal 
lines lying along the diameters of the structure and having constant angular spacing (see also fig. 7). 
There are either zero (k = 0), one (k = 1), two (k = 2), or more (k > 2) nodal diameters in bending or 
torsion modes to 73 nodal diameters because there are 146 blades on the disk. For a single 
cantilevered blade, series 1 is associated with the first natural frequency, series 2 with the second 
natural frequency, etc., where k is the number of nodal diameters. Rzadkowski et al. (2003) 
determined the influence of shaft flexibility on the natural frequencies of the shrouded bladed disks 
to four nodal diameters for two first-frequency series. Here, to four nodal diameters the number 4 is 
dependent on the geometrical parameters of the blade, disk, and shaft. Sometimes it is only one nodal 
diameter. For unshrouded cantilevered blades the natural frequencies (given in table IV) are the limit 
of the series of natural frequencies of unshrouded bladed disks. The actual frequencies are then found 
by multiplying by a constant for the configuration (not provided; contact the authors). The relative 
nondimensional amplitude values and the order of disk vibration are taken from figures of the mode 
shapes (not provided; contact the authors). 
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